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I.  INTRODUCTION 


Radio  frequency  (rf)  systems  that  utilize  the  polarization  characteristics 
of  the  target  and  environment  to  detect,  track,  identify,  etc.,  are  referred 
to  as  rf  polarimetric  systems.  These  systems  usually  combine  the  polarization 
characteristics  with  frequency  agility  for  Increased  range  resolution.  Examples 
of  such  radar  systems  are  the  >bilti-environment  Active  Radio  Frequency  Seeker 
(MARFS),  the  Advance  Indirect  Fire  System  (AIFS),  the  Helicopter  All  Weather 
Fire  Control  and  Acquisition  Radar  (HAWFCAR),  and  the  Polarimetric  Technology 
Seeker  (PTS)  as  well  as  various  other  R&D  radars  under  development  by  numerous 
independent  contractors  in  private  industry. 

The  demand  for  a  more  complete  understanding  of  the  techniques  and  pro- 
cessess  employed  in  various  programs  has  precipitated  the  development  of  the 
polarimetric  radar  simulation.  This  document  covers  the  mathematical  analysis 
required  as  background,  the  computer  simulation  model,  and  typical  results. 
Recommendations  for  future  expansions  of  this  model  are  also  addressed; 

II.  MATH  MODEL  DEVELOPMENT 

A.  Polarization  Definition 


The  concept  of  polarization  and  the  associated  conventions  are  vital  to 
the  understanding  of  the  use  of  the  polarization  scattering  matrix.  The  defi¬ 
nitions  of  polarization  have  been  traditionally  either  the  physics  or  the 
engineering  convention.  Either  convention  will  provide  the  same  general 
answer  but  with  different  notation.  Therefore,  the  convention  to  be  used 
throughout  this  analysis  is  as  stated  in  the  IEEE  STD  211-1977  "IEEE  Standard 
Definitions  of  Terms  for  Radio  Wave  Propagation". 

Linearly  Polarized  Wave  -  An  electromagnetic  wave  whose  electric  and 
magnetic  field  vectors  always  lie  along  fixed  lines  at  a  given  point.  (Page  9.) 

Left-handed  (counterclockwise)  polarized  wave  -  An  elliptically  polarized 
electromagnetic  wave  in  which  the  rotation  of  the  electric  field  vector  with 
time  is  counterclockwise  for  a  stationary  observer  looking  in  the  direction 
of  the  wave  normal. 

NOTE:  For  an  observer  looking  from  a  receiver  toward  the  apparent  source  of 
the  wave,  the  direction  of  rotation  is  reversed.  (Page  9.) 

The  definition  of  right-handed  is  found  on  page  12  and  is  the  same  as  above 
with  the  word  clockwise  used  instead  of  counterclockwise. 

B.  Plane  Waves 

For  a  plane  time  harmonic  electromagnetic  wave  traveling  in  free  space 
the  electric  field  intensity  vector  E(t),  and  the  magnetic  field  intensity 
vector  H(t)  are  always  orthogonal  to  one  another  and  have  directions  specified 
by  the  right  hand  rule  as  defined  in  the  complex  Poynting  vector  (S). 
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_  Since  E  and  H  are  always  coupled  together,  it  is  customary  to  specify  the 
E(t)  vector  only  in  describing  the  plane  wave.  The  plane  wave  can  be  spe¬ 
cified  by  its  amplitude,  frequency,  direction  of  propagation,  and  polariza¬ 
tion. 

The  vector  wave  equations  for  waves  in  free  space*  can  be  written  as 
V2E  +  k2E  -  0 
72H  +  k2H  -  0 

where  k  is  the  complex  wave  number.  The  rectangular  components  of  E  and  H 
satisfy  the  complex  scalor  wave  equation  (commonly  called  the  Helmholtz 
equation): 

V2T  +  k2¥  -0 

The  solution  to  the  Helmholtz  equation  for  one  component,  say  x,  thus  reduces 
to 

d2E*  +  k2Ex  -  0 

V 

which  is  the  one  dimensional  Helmoholtz  equation.  The  equation  has  solutions 
that  are  linear  combinations  of  e-3kz  and  e~Jkz.  We  can  choose  to  work  with 
either  of  these  solutions,  though,  in  engineering  we  generally  use  the  form 
e-jkz.  in  particular,  consider  the  solution 

Ex  -  E0e~Jkz 

This  satisfies  the  V*E-0,  and  is  therefore  a  possible  electromagnetic 
field. 

To  interpret  this  solution,  let  Eg  be  the  rms  value;  then  the  instan¬ 
taneous  field  is  found  to  be 

-yTi  E0  cos  (wt-kz) 
ly  -y/T  E0  cos  (wt-kz) 

For  conventions’  sake,  the  x  direction  will  be  the  horizontal  polarization 
and  the  y  direction  will  be  the  vertical  polarization.  In  general  the  two 
waves  need  not  have  the  same  phase.  Again,  for  convention,  it  will  be 
assummed  that  all  phase  shifts  between  the  two  waves  are  referenced  to  the 
horizontal  wave.  Therefore,  the  final  form  of  the  wave  equation  can  be  writ¬ 
ten  as 

Ex(r.t)  -  EHeJ (wt-kz )^ 

Ey(r,t)  -  EveJCvt-kz+fio)^ 

*"Time-Haroonic  Electromagnetic  Fields”  Roger  F.  Harrington,  McGraw-Hill  Book 
Co.,  1961. 
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Ej}  Is  electric  field  strength  polarized  in  the  horizontal  direction. 

Ey  is  electric  field  strength  polarized  in  the  vertical  direction, 

w  is  the  radian  frequency  of  the  transmitted  wave, 
k  is  the  complex  wave  number.* 
t  Is  time 

z  is  distance  (when  z  -  range  to  target  z  -  R) 

0O  is  the  phase  difference  between  horizontal  and  vertical  electric  field 
waves  at  the  transmitting  antenna.  (-n<^  <x) 

is  a  unit  vector  in  the  Y  direction  (vertical) 

is  a  unit  vector  in  the  X  direction  (horizontal) 

az  is  a  unit  vector  in  the  Z  direction  (range) 


C.  Special  Cases  of  Polarization 
1.  Linear  (see  Figure  1) 


Po  -  0 

VT  -  EveKwt-kz)^  or  Ey  cos  (wt-kz)ay 
HT  -  EHej(wt-kz)^  or  Eh  cos  (wt-k2)yx 
p  »  arctan  Vf/H-r  ■  arctan  Ey/E^ 


*(k«*k'-jk")  where  k'  is  the  intrinsic  phase  constant  and  k"  is  the 
intrinsic  attenuation  constant.  When  no  attenuation  is  assumed 
k  -  k’  -  2w/\. 


Figure  1  Linear  polarization 


p  equal  zero  is  referred  to  as  horizontal  polarization, 
p  equal  ninety  degrees  is  referred  to  as  vertical  polarization 
p  equal  forty  five  degrees  is  45  degree  linear  polarization. 


2.  Circular  (see  Figures  2  and  3) 

P0  -  +  90° 

By  ■  %  -  E 

Left  hand  circular  0O  -  90°  or  it/2  radians. 

The  loci  is  a  circle  of  radius  E.  The  electric  field  vector  is  constant 
in  magnitude.  When  looking  in  the  direction  of  travel  the  electric  field 
vector  rotates  counterclockwise:  when  looking  against  the  direction  of  tra¬ 
vel  the  vector  rotates  clockwise. 


I 


Figure  2.  Left  hand  circular  polarization. 


Right  hand  circular 

0O  -  -90°  or  -  n/2  radians 

r 

Ey  -  Eg  ■  E 


The  Loci  is  a  circle  of  radius  E.  However,  the  electric  field  is  rotating 
clockwise  with  time  when  viewed  in  the  direction  of  travel,  and  counterclock¬ 
wise  When  the  observation  is  made  looking  against  the  direction  of  travel. 


i  travel  of  E  is  out  of  the  page 

Figure  3.  Right  hand  circular  polarization 

Elliptical  (see  Figure  4) 

sin  P0  >  0  left  hand  elliptical 
sin  p0  <  0  right  hand  elliptical 
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angle  p  Is  the  angle  to  the  major  axis  e  and  is  dependent  upon  the  ratio 
Eu,  and  SQ. 


D.  Polarization  Notation 


As  previously  presented,  a  pure  left  hand  circular  polarized  wave 
electric  field  may  be  shown  as 

ETi  *  E[cos(wt-kz)Sx  -  sin(wt-kz)3y ) 

and  a  pure  right  hand  circular  polarization  as 

£t  =  E(cos(wt-kz)g  +  sin(wt-kz)3y] 

For  simplicity  the  time  dependency  given  originally  as  eJwt  may  be  suppressed  or 
removed  and  a  circular  wave  can  be  represented  as 

ETl  =  E[cos(-kz)ax  -  sin(-kz)ay] 

Ejr  -  E[cos(-kz)ax  +  sin(-kz)ay} 

Assuming  the  electric  field  at  the  transmitter  (z=0)  is  90°  (plus  or  minus)  out 
of  phase  in  the  H  and  V  direction,  or 

\  "  +  EeW2*y 

Ej  =  E  +  jE  for  left  hand  circular  (see  Figure  5) 

L 

Therefore 

ETr  ■  E  -  jE  for  right  hand  circular 


Y  1 

i 


Figure  5.  Left  hand  circular  wave  traveling  in  Z  direction. 


6 


E.  Scattering  Matrix 


Scattering  of  a  wave  by  objects  in  the  field  of  view  Is  modeled  by  the 
polarization  scattering  matrix  as 

1 


[E*l  -  [S][ET) 


where  E^  is  received  electric  field  vector 

^  is  transmitted  electric  field  vector 


[SJ 


c 


SHH 

SVH 


SHV 

SW 


It  can  be  shown  that  the  scattering  matrix  is  related  to  radar  cross  section 
in  the  following  manner: 


•^Ogye-J^HV 

-y/w^VV 


For  convenience,  the  —  term  is  usually  dropped  and  the  received  electric 


field  components  are  shown  to  be  related  by 


[E*]  -  [S]  [IT] 


Therefore,  for  a  monostatic  radar  the  voltage  at  the  antenna  terminals  is 
related  as 


(£«] 


l_ 

K 


where  K  is  some  factor  that  represents  the  appropriate  radar  range  scaling, 
which  for  most  calculations  is  not  considered  unless  the  absolute  received 
voltage  is  required. 

For  a  left  hand  circular  transmitted  wave  this  becomes 

EeJ(wt-2kR> 

Ee j(wt-2kR+n/2 ) 

where  R  is  now  the  one  way  range  to  the  target  from  the  radar. 
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In  short  hand  notation  this  can  be  written  as 


[Er] 


[S] 


fT- 

JET_ 


F.  Scattering  Matrix  for  Simple  Objects 


The  polarization  scattering  matrix  In  Its  most  generic  form  Is  written 
[SJ  =■ 


<J*11 

11 

S21^2, 

eJ({,12 

12 

S22ej<J>22 

where  the  subscripts  refer  to  orthogonal  components,  the  first  subscript 
being  receive,  and  the  second  transmit. 


In  the  linearly  polarized  form  this  becomes 

,S1  -  fsHHe^HH  SHVe^Hvj 

1  [_Sv„ej*VH  Swe  J 

In  the  circular  polarized  form  this  becomes 
fs]  .  [SgReJ^RR  SRLej*RL"| 

[sLReJ^LR  Sn^eJ^LLj 

where  R  refers  to  right  hand  circular,  and  L  to  left  hand  circular. 

In  this  analysis  where  a  circular  wave  (right  or  left)  Is  broken  Into 
Its  horizontal  and  vertical  components  the  linearly  polarized  scattering  matrix 
must  be  used.  However,  the  same  results  could  be  obtained  by  using  the  cir¬ 
cular  scattering  matrix  and  not  breaking  down  the  electric  field  Into  orthogo¬ 
nal  components  of  ER  and  Ey 


Consider  an  odd  bounce  reflector  (a  flat  plate)  that  totally  reflects  the 
transmitted  wave.  The  linear  scattering  matrix  elements  can  be  written  as 


S11  "  SHH»  s12  ”  SHV.  S21  "  SVH.  s22  *  SW 

The  return  from  an  Impinging  horizontal  electrical  field  will  have  the  same 
magnitude  but  the  phase  will  shift  180°.  The  same  is  true  for  an  Impinging 
vertical  electric  field.  A  horizonal  electic  field  striking  a  flat  plate  and 
being  received  in  the  vertical  direction  is  zero.  The  same  is  true  for 
transmit  vertical/receive  horizonal. 


The  scattering  matrix  for  a  flat  plate  Is  therefore 


(SJfP  ■ 


-1 

0 
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It  should  be  noted  that  the  180°  phase  shift  is  due  to  the  electromagnetic 
boundary  condition  of  zero  tangential  field  at  the  surface  of  a  perfect  con¬ 
ductor.  This  matrix  could  have  been  written  as 


fS]pp 


e-J*  0  ' 
0  e-j* 


noting  that  e  ■  -1. 

The  same  odd  bounce  reflector  in  a  circlar  scattering  matrix  would  have 
the  following  elements 

S11  “  SRR!  s12  ”  SRL;  S21  -  SLR;  S22  -  SLL 

The  rotation  of  the  return  from  a  circularly  polarized  wave  will  be  the 
reverse  of  the  rotation  of  the  transmitted  wave;  that  is,  right  hand 
transmitted  becomes  left  hand  received.  Because  there  is  only  a  cross  polari¬ 
zation  component,  the  circular  scattering  matrix  for  a  flat  plate  (odd)  is 


(SIfp 


P  1 

Li  oj 


The  equivalence  of  linear  polarization  and  circular  polarization  can  best 
be  seen  from  examples  of  both  worked  in  parallel.  Assume  a  left  hand  circular 
transmit  into  a  flat  plate. 

The  circular  (left  hand)  transmit  is  written  in  a  linear  system  as 
-  Eyx  +  jE^y 

A  circular  transmit  system  Is  written  In  circular  notation  as 
T  T 

eT5t  -  Er 

where  R  is  a  unit  vector  rotating  in  the  right  hand  direction 
L  is  a  unit  vector  rotating  in  the  left  hand  direction 


LINEAR  (FLAT  PLATE) 


CIRCULAR  (FLAT  PLATE) 


Noting  that  the  direction 
of  travel  has  reversed 
the  received  wave  is  of 
the  form 


which  is  a  right  hand 
circular  wave  traveling  in 
the  -z  direction. 


In  the  circular  form 
the  received  wave  is 
of  the  form 

E  »  e£r  +  0 

which  is  a  right 
hand  circular  wave. 


Figure  6.  Scattering  characteristics. 


The  analytical  relationships  developed  for  the  simulation  are  based  upon 
the  linear  scattering  characteristics  of  a  few  simple  shapes,  classified  to 
some  degree  by  the  number  of  reflecting  surfaces  encountered. 

1.  Odd  bounce  scattering  matrix  (flat  plate,  trihedral  corner  reflector) 


for  linear  polarization  (see  Figure  6) 


2.  Even  bounce  scattering  matrix 
cos2A  sin2Q 

sin20  -cos20 

where  9  is  the  angle  of  rotation  of  th< 
(See  Figure  7). 


(diplane)  for  linear  polarization 


diplane  relative  to  the  horizontal. 
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w  ^ 


Figure  7.  Diplane  rotation  angle. 

3.  Dipole  matrix  for  linear  polarization  (Figure  8). 

"-cos  20  -cosSsinS 

-cos0sin0  -sin20 


( 


Figure  8.  Dipole  rotation  angle. 

It  is  assumed  that  superposi tion  holds  such  that  a  complex  target  may  be 
modeled  by  an  ensemble  of  these  even  and  odd  bounce  targets  or  scatterers  with 
the  inclusion  of  their  respective  ranges. 

G.  Radar  Range  Scaling 

The  basic  radar  equation  to  determine  the  power  received  at  the  radar  is 

„  PfG2X2o 

Pr  “ 

i 

where  Pr  is  power  received  in  watts 

Pt  is  peak  power  transmitted  in  watts 
G  is  antenna  gain  (unitless) 

X  is  wavelength  in  meters 
o  is  radar  cross-section  in  meters  squared 
R  is  range  to  target  in  meters 
Ls  is  system  loss  (unitless) 

Because  this  analysis  is  performed  in  the  voltage  domain  the  standard 
radar  equation  must  be  modified  to  be  expressed  in  terms  of  voltage. 
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where  Vpeak  Is  peak  voltage  received 

Z  Is  Impedance  (assumed  50  ohms) 

Therefore,  the  peak  voltage  (Vpeak)  Is 


Vpeak  -  V 2*Z*Pr 


By  removing  the  radar  cross-section  from  Pr,  Pr  becomes  a  constant  radar 
scalar  which  when  used  with  the  field  voltage  obtained  from  the  scattering 
matrix  defines  the  peak  received  voltage. 


Vpeak 


'v! 


2ZPtG2X2 


(An)1]^, 


jr 


H.  Noise  Generation 

If  the  radar  were  operated  In  a  perfectly  noise  free  environment  so  that 
no  external  noise  sources  accompanied  the  desired  signal,  there  would  still 
exist  an  unavoidable  component  of  noise  generated  by  the  thermal  motion  of  the 
conduction  electrons  In  the  receiver  input  stages.  This  Is  called  thermal 
noise  and  Is  directly  proportional  to  the  temperature  of  the  ohmic  portions  of 
the  circuit,  and  the  receiver  bandwidth.  The  available  thermal-noise  power 
generated  by  a  receiver  of  bandwidth  Bjj  (In  Hz)  at  temperature  T  (degrees 
Kelvin)  Is  equal  to: 

average  available  power  *  KTBn 
where  K  Is  Boltzmann's  constant  (1.38  x  lO-2^  joule/deg) 

No  matter  whether  the  noise  is  generated  by  a  thermal  mechanism  or  by 
some  other  mechanism,  the  total  noise  at  the  output  of  the  receiver  may  be 
considered  to  be  equal  to  the  thermal-noise  power  obtained  from  an  ideal 
receiver  multiplied  by  a  factor  called  noise  figure  (NF).  The  noise  figure 
(NF)  of  a  receiver  Is  defined  by  the  equation: 

NF  »  noise  out  of  practical  receiver 

noise  out  of  ideal  receiver  at  Std  Temp  (T0) 

The  standard  temperature  is  taken  to  be  290°  K. 

Therefore, 


average  available  power  =*  KT0BNF 

Assuming  this  to  be  the  available  average  power  at  the  input  stages  of  a 
radar,  the  ohmic  load  is  assumed  to  be  matched  as  In  the  simple  circuit 
diagram  in  Figure  9. 


Figure  9.  Equivalent  noise  circuit. 


NOTE:  Load  resistance  (Rl)  Is  matched  to  source  resistance  (Rs). 

Therefore,  the  RMS  voltage  available  at  the  source  can  be  calculated  as 


V8  -  2VL 

Vl  ~^KT0BNFRl  and 


!\AtobnfrL 


The  noise  entering  the  IF  amplifier  Is  assumed  to  be  Gaussian,  with  a 
probability-density  function  given  by 


p(v)dv  - 


*{£) 


dv 


where  p(v)  dv  Is  the  probability  of  finding  the  noise  voltage  between  the 
value  of  v  and  v  +  dv,  T0  Is  the  variance,  or  mean-square  value  of  the  noise 
voltage.  The  mean  value  of  v  is  taken  to  be  zero. 

2 

Therefore,  the  mean-square  value  la  taken  to  be  or  KT0BNFR^  and  the 
standard  deviation  by  definition  Is 

SD  «y'KT0BNFRL 
I.  Antenna  Isolaton 


When  two  antennas  (or  elements)  are  widely  separated  the  energy  coupled 
between  them  is  small,  and  the  Influence  of  the  receiving  antenna  on  the 
current  excitation  and  pattern  of  the  transmitting  antenna  is  negligible.  Ai 
the  antennas  (or  elements)  are  brought  closer  together  the  coupling  between 
them  increases. 

Isolation  of  a  polarimetrlc  antenna  Is  represented  as  two  antennas  that 
cross  couple  energy  during  transmit  and  receive.  Considering  only  the 
transmit  cycle  the  coupling  can  be  represented  as  In  Figure  10. 


Figure  10.  Transmit  isolation. 


In  Figure  10  Ey  and  Eg  are  input  to  the  antenna,  and  V-p  and  H-p  are  antenna 
outputs. 


r2  „  10-ISOL/20 


Rp  yi-ioisM7io 


ISOL  is  the  antenna  isoluatlon  in  dB  one  way,  always  positive. 

Assuming  reciprocity,  the  isolation  upon  receive  is 
eVR  “  vRr1+hRr2 
eHR  “  hRr1+vRr2 

where  Eyg  and  Egg  are  the  input  signals  to  the  receiver  and  Hg  and  Vg  are 
the  inputs  at  the  antenna  plane. 

Phase  stability  is  assumed  across  the  antenna  plane. 

J.  Frequency  Agility  and  Intra-Range  Resolution 

Range  resolution  is  usually  defined  as  the  distance  at  which  two  targets 
can  be  resolved  in  range.  In  the  conventional  radar  this  is  defined  by  the 
pulse  width  of  the  transmitted  wave  as  AR  »(ct/2) 

where  AR  is  range  resolution  (m) 

C  is  velocity  of  light  (m/sec) 
x  is  radar  pulsewidth  (sec) 

Considering  the  radar  to  have  a  match  receiver  1,  equal  to  one  over  receiver 
bandwidth,  AR  becomes 

“'ft 

where  B  is  bandwidth  in  Hertz. 

Either  of  the  two  equations  can  be  used  to  calculate  the  range  resolu¬ 
tion  of  a  radar.  However,  the  latter  equation  is  the  more  general  form  and 
can  be  ulltized  in  calculating  range  resolution  in  conventional  radar,  pulse 
compression  radar,  and  frequency  agile  radar,  as  well  as  in  hybrids  of  these 
such  as  the  pulse  compression  frequency  agile  radar. 

Ruttenberg  showed  in  1967  a  method  that  increased  range  resolution  with 
a  non-coherent  source.  This  involved  a  frequency  agile  scheme  that  summed  the 
pulses  after  they  were  received  (coherent  on  receive)  and  delayed  by  1/PRF. 

Since  then  the  use  of  a  fully  coherent  radar  utilizing  frequency  agility, 
pulse  to  pulse,  and  the  Digital  Fourier  Transform,  has  demonstrated  a  range 
resolution  technique  that  does  not  require  delay  lines  as  did  Ruttenberg 's 
technique.  The  coherent  pulses  are  fed  to  a  DFT  (usually  the  same  size  as  the 
number  of  frequency  shifts)  and  frequency  is  transformed  into  time  (via  the 
DFT)  with  intra-range  resolution  of  the  system  following  the  same  range  reso¬ 
lution  equation. 


:  -  A. 


AR 


C 

2B 


where  B  Is  now  the  frequency  agile  bandwidth. 

G jessing.  In  his  book  "Adaptive  Radar  In  Remote  Sensing"  shows  that  the 
amplitude  spectrum  of  the  scattered  field  is  the  Fourier  transform  of  the 
delay  function  f(t).  Thus,  If  the  target  is  at  some  distance  d,  the  delay 
function  will  oscillate  with  a  period  c/2d.  Therefore,  by  the  use  of  a 
multi  frequency  radar  system,  the  resolution  of  the  radar  can  be  increased  as 
the  bandwidth  of  the  agile  radar  increases.  The  complex  Fourier  transform 
will  provide  the  true  reference,  while  the  amplitude  only  Fourier  transform 
will  provide  the  relative  distances  between  the  resolvable  elements. 

III.  RADAR  SIMULATION 

The  functional  diagram  of  a  polarimetric  radar  is  shown  in  Figure  11. 
Functionally  the  model  is  a  frequency  agile  coherent  radar  model.  If  a  non¬ 
coherent  radar  model  is  desired  the  signal  processor  section  can  be  modified. 
The  frequency  agile  waveform  selects  the  transmitter  and  coherent  local 
oscillator  frequency.  The  transmitter  energy  is  split  (coupled)  to  the  dual 
polarized  antenna  with  a  +90°  phase  shifter  in  the  vertical  channel,  resulting 
in  either  right  hand  or  left  hand  circular  polarization.  Zeroing  one  channel 
or  the  other  results  in  horizontal  or  vertical  polarization.  Removing  the 
phase  shifter  and  adjusting  the  splitter  results  in  slant  polarization  of  any 
desired  angle.  If  the  cross  coupling  in  the  antenna  section  is  large  enough 
the  result  is  an  elliptical  wave. 

Energy  reflected  from  the  target  area  is  received  in  both  the  horizontal 
and  vertical  antennas  with  cross  coupling,  and  passed  to  the  coherent  Inphase 
(I)  and  Quadrature  (Q)  detectors,  resulting  in  IF  detected  signals  of  horizon¬ 
tal  I  and  Q,  and  vertical  I  and  Q. 

This  radar  simulation  configuration  allows  maximum  versatility  by  pro¬ 
viding  for  circular,  elliptical,  and  linear  polarization  transmission. 
Receiving  horizontal  and  vertical  with  antenna  cross  coupling  allows  the 
signals  of  circular  and  elliptical,  and  horizonal  and  vertical,  either 
coherent  or  non-coherent.  The  configuration  shown  in  Figure  11  is  not 
Intended  to  imply  prefered  hardware  configuration,  but  rather  to  depict  a 
radar  simulator  which  can  be  used  to  simulate  a  large  number  of  pulsed  polari¬ 
metric  radars  in  order  to  evaluate  proposed  radars  and  signal  processors. 

A.  Signal  Processing 

1.  General 

Outputs  from  the  radar  simulation  (HI,  HQ  and  VI,  VQ)  are  input  to  the 
signal  processing  software  where  they  are  combined  to  form  various  types  of 
received  signals  and  the  respective  inverse  Fourier  transforms.  The  radar 
signals  available  from  the  the  signal  processor  as  plots  are: 
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Figure  11.  Radar  configuration. 


a.  Peak  Horizontal  voltage 

b.  Peak  Vertical  voltage 

c.  Phase  between  Horizontal  and  Vertical 

d.  Peak  Left  Hand  Circular  (LHC)  voltage 

e.  Peak  Right  Hand  Circular  (RHC)  voltage 

f.  Phase  between  LHC  and  RHC 

g.  Scatterer  Locations  (Inverse  FFT's  of  coherent  signals) 

h.  Scatterer  Seperatlon  (Inverse  FFT's  of  non  coherent  signals) 

1.  Phase  plots  of  FFT's 

2.  Linear  Polarlzaton 

Functionally  the  coherent  horizontal  and  vertical  signals  are  processed 
as  shown  In  Figure  12.  Hie  resulting  lnphase  and  quadrature  signals  are  then 
loaded  Into  a  complex  array,  and  an  Inverse  FFT  is  performed.  The  resulting 
lines  represent  the  location  of  the  scatterers  relative  to  the  leading  edge  of 
the  radar  range  gate. 


Quadrature  Detector  (Double  Balanced  Mixer) 

Horizontal  Received  Signal  »  HI+jHQ 
Vertical  Received  Signal  -  VI+jVQ 


Figure  12.  Linear  coherent  detection  block  diagram. 

Loading  the  real  portion  of  the  complex  array  with  the  amplitude  magni¬ 
tudes  only,  and  performing  an  inverse  FFT,  results  in  lines  that  represent  the 
separation  of  scatterers  relative  to  each  other.  There  is  one  line  (neglecting 
sldelobes)  for  each  combination  of  pairs  of  scatterers. 

N 

No.  lines  -  T.  (»-i) 
i-1 

where  N  is  the  number  of  reflectors. 
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3.  Circular  Polarization 


Figure  13  is  a  functional  block  diagram  of  a  linear  to  circular  tranfor- 
mation.  Inputting  horizontal  and  vertical  lnphase  and  quadrature  results 
in  left  and  right  hand  circular  polarized  signals.  This  can  be  repre¬ 
sented  In  matrix  notation  as: 


H-vf-  [!'!]  [I] 


Loading  the  resulting  RHC  or  LHC  into  a  complex  array,  and  performing  an 
Inverse  FFT,  results  In  lines  that  represent  the  location  of  scatterers  rela¬ 
tive  to  the  leading  edge  of  the  radar  range  gate.  ^ 

Loading  the  real  portion  of  the  complex  array  with  the  amplitude  magnitudes 
only,  and  performing  an  inverse  FFT,  results  in  lines  that  represent  the 
seperatlon  of  scatterers  relative  to  each  other. 

IV.  SIMULATION  ULITIZATION 

A.  General  Simulation  Outputs 

Figures  14  through  35  are  the  22  output  plots  of  the  simulation.  Four 
received  amplitudes  (Horizontal,  Vertical,  LHC  and  RHC)  are  plotted  as  a  func¬ 
tion  of  the  transmit  frequency.  The  header  data  presented  at  the  top  of  each 
plot  define  the  radar  operating  parameters  used  in  the  calculations.  All 
plot  headers  for  each  input  parameter  set  are  identical.  The  phase  angle  bet¬ 
ween  received  Horizontal  and  Vertical  or  RHC  and  LHC  is  plotted  as  a  function 
of  transmit  frequency  and  is  the  angle  Beta  discussed  in  paragraph  2.2.  The 
remaining  16  plots  are  inverse  FFT's,  amplitude  and  phase,  plotted  as  a  func¬ 
tion  of  intra-range  gate  resolution.  The  FFT  amplitude  plots  labeled  I&Q  pro¬ 
vide  scatterer  location  from  the  leading  edge  of  the  range  cell  while  the 
plots  labeled  peak  amplitude  provide  scatterer  separation  depending  on  whether 
the  amplitude  data  were  loaded  into  the  FFT  as  complex  ISQ  or  as  real  ampli¬ 
tude  only.  FFT*s  were  loaded  in  ascending  order  with  the  received  signal  from 
the  lowest  transmit  frequency  in  location  one. 
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B.  Antenna  Isolation 


Utilizing  the  simulator  program,  and  varying  the  amount  of  one  way 
antenna  polarization  isolation,  can  reveal  the  effects  of  Isolation  on  the 
polarlmetric  outputs  of  a  system.  This  Is  exampled  by  Figures  28  through  32 
which  show  the  LHC  and  RHC  scatterer  locations  for  a  four  target  array  with  30 
dB  one  way  polarization  Isolation.  Figures  36  and  37  present  the  same  con¬ 
ditions  with  10  dB  Isolation  for  comparison.  Comparing  these  plots  one  can 
observe  the  cross  coupling  from  one  channel  to  the  other. 

While  this  example  is  presented  for  LHC  and  RHC  output  it  is  obvious  that 
the  other  outputs  of  the  simulation  may  also  be  examined.  Antenna  Isolation 
effects  on  horizontal  or  vertical  outputs.  In  either  the  coherent  or  non¬ 
coherent  mode,  as  well  as  other  combinations  of  transmitter  and  receiver  polari¬ 
zation  configurations,  can  be  examined. 

C.  Signal  to  Noise  Ratio 

System  noise  effects  on  polarlmetric  outputs  can  be  examined  in  two  ways: 
first,  by  increasing  target  range  (reducing  signal  strength);  second,  by 
increasing  the  receiver  noise  figure  (Increase  system  noise).  Examples  of 
these  are  given  in  Figures  17  and  21  (horizontal  and  vertical  scatterer  loca¬ 
tions  for  a  greater  than  30  dB  signal  to  noise  ratio),  and  Figures  38  and  39 
(for  a  signal  to  noise  ratio  of  8  dB). 

Inputting  a  clutter  model  and  varying  the  system  noise  will  allow  examina¬ 
tion  of  the  effects  of  the  clutter  to  noise  problem  on  signal  processing. 

System  noise  can  be  increased  by  elevating  the  receiver  noise  figure. 

D.  Signal  to  Clutter 

The  utilization  of  the  simulator  program  to  explore  the  effects  of  clutter 
on  signal  detection  will  be  highly  dependent  upon  the  target  and  clutter 
models  used.  A  model  for  clutter  in  polarlmetric  form  that  has  been  truly 
verified  has  yet  to  be  developed.  Therefore,  in  order  to  demonstrate  the  use  of 
the  program  the  following  example  will  be  used:  a  contrived  target  model  of  one 
and  one  half  meters  radar  length,  made  up  of  five  reflectors  randomly  spaced, 
and  having  a  radar  cross  section  of  five  square  meters  each  (Figures  40  through 
43);  clutter  made  up  of  fifty  randomly  selected  location,  orientation,  and  type 
spaced  reflectors  of  0.1  square  meters  each.  This  example  has  a  total  signal  to 
clutter  ratio  of  23/5,  or  7  dB,  and  Is  shown  in  the  horizontal  and  vertical 
location  plots  in  Figures  44  and  47.  Figures  48  through  51  show  the  same  con¬ 
figuration  with  the  clutter  cross  section  Increased  to  one  square  meter  per 
reflector.  This  represents  a  signal  to  clutter  ratio  of  -3.0  dB  (25/50). 

V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  digital  simulation  has  been  developed  to  investigate  various  aspects  of 
a  frequency  agile,  polarlmetric  pulsed  radar  system.  While  the  simulation  Is 
not  all  inclusive  and  will  undoubtedly  be  refined  and  updated  for  years  to 
come,  it  is  a  useful  tool  for  evaluating  both  hardware  and  software  effects  on 
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the  next  generation  of  redars.  The  simulation  validation  was  performed  by 
comparison  with  an  operational  radar.  The  validation  has  been  excluded  as  the 
data  were  acquired  from  a  contractor’s  IR&D  radar.  Any  government  agency 
desiring  more  information  relative  to  the  validation  should  contact  the 
authors  at  AV  746-4061. 

Major  limiting  factors  to  simulation  results  are  the  target  and  clutter 
scattering  models  which  remain  basically  undefined  at  this  time.  It  is  recom¬ 
mended  that  all  models  and  data  in  the  future  be  done  in  scattering  matrix 
format  so  that  the  entire  radar  signature  will  be  available  for  future  radar 
hardware  and  simulator  designers.  Without  such  data  and  validated  models  the 
radar  system  analyst  and  designer  will  continue  to  suffer  from  the  so  called 
"Sedenqulst  Effect";  that  is,  put  two  radar  engineers  in  a  room  and  say  the 
word  "clutter";  return  years  later  and  they  will  still  be  trying  to  define  and 
agree  as  to  what  clutter  is. 

Future  plans  for  the  polarlmetric  radar  simulator  include  the  addition  of 
doppler,  tracking  errors,  jamming,  attenuation  back  scatter  due  to  weather, 
and  realistic  model  development. 

This  simulation  does  not  include  cross  range  positional  effects  of  scat- 
terers.  All  cross  range  scatterer  positions  within  the  antenna  beamwidth  are 
collapsed  to  a  single  radial  range  bin.  The  Inclusion  of  doppler  will  pro¬ 
vide  the  second  "cross-range"  dimension  for  two  dimensional  analysis. 
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FREQUENCY  ^GIGAHERTZ* 


GATE  DISTANCE  (METERS) 


INTRA-RANGE-GATE  DISTANCE  (METERS) 

Figure  18.  Inverse  FFT  phase  angle  for  horizontal  I&Q  at  30  dB  antenna  isolation 


<s> 


CL  LU  <E  W  IOQJwN  <CCOL 


MZ3  U.  Ll.  f— 


26 


INTRA-RANGE-GATE  DISTANCE  (NETERS) 

Figure  21.  Inverse  FFT  of  vertical  I&Q  at  30  dB  antenna  isolation 
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INTRA-RANGE-GATE  SEPERATION  (PIETERS) 

Inverse  FFT  of  peak  vertical  amplitude  at  30  dB  antenna  isolation 


INTRA-RANGE-GATE  SEPERATION  (METERS) 

Figure  24.  Inverse  FFT  phase  angle  of  peak  vertical  amplitude  at  30  dB  antenna  isolation 
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FREQUENCY  XGIGAHERTZ* 

Figure  25.  Peak  RHC  amplitude  vs.  frequency  at  30  dB  antenna  isolation 
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itude  vs.  frequency  at  30  dB  antenna  isolation 


FREQUENCY  XGIGAHERTZ* 

Figure  27.  Phase  angle  between  RHC  and  LHC  at  30  dB  antenna  isolation 


INTRA-RANGE-GATE  DISTANCE  (METERS) 

Figure  28.  Inverse  FFT  of  RHC  I&Q  at  30  dB  antenna  isolation 


INTRA-RANGE-GATE  DISTANCE  (METERS) 


INTRA-RANGE-GATE  SEPERATION  (NETERS) 

Figure  30.  Inverse  FFT  oi  peak  RHC  amplitude  at  30  dB  antenna  isolation 
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INTRA-RANGE-GATE  SEPERATION  (PIETERS) 

Figure  35.  inverse  FFT  phase  angle  of  peak  LHC  amplitude  at  30  dB  antenna  isolation 


INTRA-RANGE-GATE  DISTANCE  < METERS ) 

Figure  36.  Inverse  FKT  of  RHC  I&1J  at  10  dB  antenna  isolation. 


INTRA-RANGE-GATE  DISTANCE  ( METERS) 

Figure  37.  Inverse  KFT  of  LHC  I&Q  at  10  dB  antenna  Isolation 


Inverse  FFT  of  horizontal  I&Q 
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^■rtical  voltage  for  target  model  only 


INTRA-RANGE-GATE  DISTANCE  ( PIETERS ) 

Figure  43.  Inverse  FKT  of  vertical  I&Q  for  target  model  only 
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Figure  45.  Peak  vertical  voltage  for  signal  to  clutter  ratio  of  +7  dB. 
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figure  46.  Inverse  FFT  of  horizontal  I&Q  tor  signal  to  clutter  ratio  of  +7  dB. 
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APPENDIX 


Simulation  Flow  Charts  and  Program  Listing 

The  following  flow  charts  were  developed  as  an  aid  In  following  the  mathe¬ 
matical  development  of  signals.  The  program  listing  of  subroutines  Is  short 
enough  to  provide  an  easily  followed  path  without  flowcharts.  This  simulation 
has  been  developed  and  run  on  a  Digital  Equipment  Corporation  (DEC)  28K  word 
PDP-11/10  computer  running  DEC'S  RT-11  operating  system.  The  plots  were  per¬ 
formed  utilizing  a  Tektronix  terminal  4014  driven  by  in-house  developed 
plotting  software.  The  plotting  subroutines  are  described  by  function  only 
without  software  listing  Included  In  the  program  printout.  This  will  provide 
a  programming  guide  for  tailoring  plots  to  other  systems. 

DEC's  RT-11  Subroutines 


Call  Assign  -  Attaches  a  disk  file  for  reading  or  writing  and  assigns  a 
logical  unit  number. 

Call  Close  -  Closes  an  attached  disk  file. 

In-House  Computer  Subroutines 

Call  LAND  -  Performs  logical  bit  anding  of  the  two  arguments. 

Call  SWR  -  Read  computer  switch  register.  Used  to  control  line  printer 
and  hard  copy  functions. 

Call  M.0GN  -  Perform  forward  or  inverse  in  place  FFT  of  complex  array. 
In-House  Plotting  subroutines 

Call  PLOT  -  Erase  4014  screen. 

Call  V14CSZ  -  Select  size  of  Alphanumeric  characters  typed  on  4014. 

Call  AXES  -  Draw  plotting  axes  by  screen  position  and  tic-marks  controlled 
by  user  units  and  store  parameters  for  user  units  plotting  by  call  LINE. 

Call  LINE  -  Draw  graphic  line  on  4014  between  two  points  described  by  user 
units . 

Call  HRDCPY  -  Cause  hard  copy  of  4014  screen  to  be  produced. 

Call  LABEL  -  Provide  X  and  Y  axes  labels  to  be  centered  and  typed  on  4014. 


Call  STALL  -  Cause  computer  to  wait  momentarily  for  HRDCPY  to  be  executed. 


INITIALIZE  CONSTANTS 

t  ,  2  it  ,(2u  )2,(4tt  )3  (1 fjt),  C,  COR,  UVOLTS 


4 

r 

INPUT 

AISOL 

— 

ANTENNA  ISOLATION 

INCREMENT  IRUN 
IRUN=IRUN+1 


CALL  IXMIT  SUBROUTINE 
GET  E1H,  E1V,  PH I  TV,  NXMIT 


R1=SQRT  ( 1 . -10.**( -AISOL/IO. ) ) 

Rl=10. **('*AI SOL/20 } CALCULATE  ANTENNA 
CROSS-COUPLING  RATIOS.  Rl-COUPLED  ENERGY 
Rl-REMAINING  ENERGY 


ANGHT-0+10 
ANGVT=0+, jPHITV 


EHTP=ElH*e ( ANGHT) 
EVTP=ElV*e(ANGVT) 


MULTIPLY  BY  CROSS  COUPLING  FOR  TRANSMIT 

EHT=EHTP*Rl+EVTP*R2 

EVT=EVTP*R1+EHTP*R2 


PHASE  SHIFT  IN  RCVR 
PHIRH=0. 

PHIRV=0. 


1 

—  ■ 

iNPUTf 

AGILE 

NSTEP 

:REQ 

STEPS 

* 

TNPUTF 

SIZE 

NFFT 

‘FT 

“Input- 

DATA  F 
IFILE 

[LE  NAME 

C 

■) 

6 


2 


EAD  DATA  FILE 
NOISE-ON  OR  OFF  (1  or  0) 

RIFBW-RECEIVER  IF  BANDWIDTH 
RNFDB-RCVR  NOISE  FIGURE  IN  DB 
CF-XMIT  CENTER  FREQ  (HERTZ) 

FBW-FREQ  AGILE  BANDWIDTH  (HERTZ) 

PRF-PULSE  REP.  FREQ  (HERTZ) 

PTPWR-XMIT  POWER  (WATTS) 

CR-COMPRESSION  RATIO 
GAINA-ANTENNA  GAIN  (dB) 

RANGE-RANGE  TO  LEADING  EDGE  OF  RANGE  CELL  (METERS) 
DBLOSS-SYSTEM  LOSSES  (dB) 

ASCALE-AMPLITUDE  PLOTTING  AMPLITUDE  SCALE  (  VOLTS) 
NSCAT-NUMBER  OF  SCATTERERS 


CALCULATE 

DF =FBW/FLOAT ( NSTEP-1 )-FREQ  STEP  SIZE 

ANTG=10.**  (GAINA/10. )-ANTENNA  GAIN 

ANTG2=ANTG**2-ANTENNA  GAIN  SQUARED 

RANGE4=RANGE**4. -RANGE  TO  THE  FOURTH  POWER 

SL0SS=10.**(DBL0SS/10.)SY$TEM  LOSSES  IN  EITHER  H  OR  V  CHANNEL 

VARI«50.*1.38E-23*290.*RIFBW*RNF 


-READ  SCATTERERS'  DATA 
TYPE,  SIZE,  ROTATION  ANGLE,  DISTANCE 


CONVERT  ROTATION  ANGLES  FROM  DEGREES  TO  RADANS 


SET  NOISE  FREQ  ACCUMULATORS  TO 

ASHRH=0 

ASNRV=C 


■ 


FIRST  FREQ  STEP  TO  1 
IF  STFQ= 1 


GENERATE  FREQUENCY  F 
FROM  FUNCTION  FREQ  (IF) 


CALCULATE  WAVELENGTH 
X  =  C/F 


CALCULATE  PHASE  SHIFT  DUE  TO  FREQUENCE 
PHIF=4  tt  /  A 


CLEAR  COMPLEX  ACCUMULATORS 

F.HACC=0+j0 

EVACC=0+jO 


63 


>-  .  -vrM 


CALCULATE  TOTAL  PHASE  SHIFT  OUE  TO  FREO  &  DISTANCE 
PHIFD=PHIF*SCATER  (1,4) 


l 


f 

CALL  A 
BY  SCA1 

SCATTERING  MATRIX  SUBROUTINE 
fTERER  TYPE 

*  ■  -  - 

> 

_  5 

* 

CALCULATE  TOTAL  PHASE  SHIFT  AS  COMPLEX  EXPONENT 
4NGH=0+j  (PHIRH+PHIFO) 

AN6V=0+1  (PHIRV+PHIFD) 

< 

f 

CALCULA 

ETHOR=E 

ETVERT= 

TE  PHASF  SHIFTED  RE-RADIATE  H&V  WAVES 

HT*e_(ANGH) 

EVT *e~<4NGV) 

_ *- 

CALCULATE  AMPLITUOE  ATTENUATION  FROM  FUNCTION  RSCALE 
SCALE=RSCALE(LAMDA) 


? 


CALCULATE  RE-RADIATED  SCATTERER  ENERGY  RECEIVED  AT  ANTENNA 
EH=SCALE  (Sn*ETHOR+S12*ETVERT) 

EV=SCALE  (S21*ETHOR+S22*ETVERT) 


SAVE  EH 
EHTP=EH 


— 

CALCUL 

EH=EH* 

EV=EV* 

! - 

ATE  ANTENNA  CROSS  COUPLING  EFFECT 
R1+EV*R2 

R1+EHTP*R2 

>  ... 

ACCUMU 

EHACC= 

EVACC= 

LATE  SCATTERER  EFFECT  IN  H&V 

EHACC+EH 

EVACC+EV 

NOISE  ON? 


NO 


GET  GAUSSIAN  NOISE  AMPLITUDE 
CALL  FUNCTION  VNOISE 
VNl=VNOISE(SD) 

VN2=VN0ISE(SD) 


GET  I&Q  OF  HORIZ  CHNL 

HI=REAL(EHACC) 

HQ=AIMAG(EHACC) 


GET  I&Q  OF  VERT  CHNL 
VI=REAL(EVACC) 

VQSA IMAG ( EVACC ) 


CALCULATE  HORIZ  RCV  SIGNAL  PHASE  ANGLE 
HZANG=ATAN2(HQ,  HI) 


CALCULATE  VERT  RCV  SIGNAL  PHASE  ANGLE 
VTANG-ATAN2(VQ,  VI) 


CALCULATE  PHASE  ANGLE  BETWEEN  H&V  REFERENCE 

TO  HORIZONTAL 

BETA-VTANG-HZANG 


NORMALIZE  BETA  TO  +  180  DEGREES 


STORE  PHASE  ANGLE  BETA 
BETSAV(IF)*BETA 


9 


69 


•V**' 


l-. 


ALL  N 
FREQUENCE 
TRANSMUTE 


NO 


D 


CALCULATE  INTEGER  POWER  OF  2  FOR  FFT  SIZE 
N-AL0GI0(FL0AT(NFFT)/AL0G10(2. )+0.5 


I 


INITIALIZE  SIGNAL  TO  NOISE  VALUES 

SNRH=99.99 

SNRV=99.99 


CALCULATE  AVERAGE  SIfiNLE  PULSE  SIGNAL  TO  NOISE  RATIO  IN  dB 
SNRH-20 . *AL0G10( 1 / VX*ASNRH/NSTEP/SD 
SNRV-20 .  *AL0G10  ( 1  f  y7*ASNRH/NSTEP  /SO 


m 


SET  UP  AUTOMATIC  HARDCOPY 
AND 

DATA  PRINTING  IF  DESIRED 


o* 


HORIZONTAL  -  1ST  ITERATION 
RHC  -  2ND  ITERATION 


VERTICAL  -  1ST  ITERATION 
LHC  -  2ND  ITERATION 


LINEAR  PHASE  ANGLE  -  1ST  ITERATION 
CIRCULAR  PHASE  ANGLE  -  2ND  ITERATION 


PLOT  PHASE 
VS  FREQ 
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PAbE  001 


OUu  1 

c 

c 

c 

C 

C 

C 

C 

C 

0002 
U0o3 
UUU  4 
0001 
OOUb 
0007 
uOuo 
oOuS 
0010 
00  11 
0ol2 
00  1  3 
0<i  l  4 
OUll 
u016 
0017 
0018 
ooio 


oo«:0 


UUc  1 

c 

c 

o0e2 
002  J 

O0<4 
0021 
0«2b 
u02  7 
oocd 

O02O 

oo  jo 
00  i  1 

c 

00  32 
0033 
0034 
0031 
003b 
0037 
0038 


PPOGhAm  SlhPTS 

PWubh  Ah  70  SIMULATE  1  nt  hF  GUIDANCE  TECHNOLOGY  •  S  POLAPImEThIC 
IEChnulOGY  SEEkEp  (PTS). 

akihe.*  or:  p.  *.  stuENuuisr  and  «.  F.  pussell  4-peo-82 

LAlESf  UPDATES  28-0CT-82 

InFEuEk  nXAXIS12o)  .NYAXISF20) 

1NFEGEP  IFILE(dJ,NXMIU2) 

PEAL  L3E I  A , L AMuA 
COMPLEX  V  ALOE  125b) 

COMPLEX  AVbVAL 

DImEnSIOn  SFPEo(21bJ 

DIME. '.31  On  but  r  AH  1 21b)  ,  SUE  T  A  V  1256) 

DIMENSION  SC A  TEH (100,4) , AHSAV (2161 ,AVS AO (25b 1 ,bETSAV(21b) 

DIr.E.\SlO.»  HwoAv  (21b)  ,  rt  I SA  V  (21b)  ,  VuSAV  (2lb)  ,  0  1  S  A  V  (21b) 
ol MENS  I UN  CUET«Sl21b) 

OImEnSIUn  A ( 20 1 

COMPLEX  PMCSA01236) , L nCSa V 125b ) 

COmPlex  34ATitXi2«2),cM»Ev,ETnON,ETvEPF,ANGH,ANGV,£HTP,tVTP,EHT,£vF 
COMPLEX  AnGMT.A.MGVT 
COMPLEX  EMACC.EVACC, ONH, VNV 

common  /NOKKr /1FSTE0, IUP,LSTEP,N8TEP»0F,CF,Ebh 

COMMON  /nKSC I /  SCA  T EH  v  SMA  TP X 

COMMON  /HE AO/ A I SUL,. NSC A T,GA I NA, NOISE. 

I  P  AnGE , UBLOSS, nXMI I , 1F1LE.SU, 816, 

1  SNPn, SnHV, SnRnI , SNPMJ, SNKO I , SNkVQ, S.«R 
COMMON  /aUn*/muSAV,HISAV,vuSAV, VISAV,SM»tu,SbETAH,5bt Fay.BETSAV, 

1  PHCS A  V .LnCSAO.CSt TaS, AhSAV , AvSAV, 

1  VAI.UE 

COMMON  /Sib  ,AL/P  TPaP,  KAnGE4,Cp ,  AnTG2,  SLUSS,  P14C 

INI HaLI/ATIOn  OALUtS 
V 1USn2= 1 . /SOP T  l2 ,  ) 

uVOLFS=l.E-b  JMICRO-VoLTS  SCALE* 

C  =  2.'H  743ES 
1CPY=0 
iPo.»s  1 
P 1=3,14119 
Pl2=2, *P  1 
PI4=2.«P12 
PI4C=PI4»*3, 

COp  =  p I / 1  so .  JCunvEPT  DEGREES  10  RADIANS 

call  plotio) 
call  V14CS/U) 
type  bool 

ACCEPT  bO  04 , A [ oOL 

CALL  XmII (EIh.eIV.PhI I V.NXMl 1 » 

P  l=SuP  f  1 1  .-l  o.«*  1-AISoL/  10,  )  J  !PE4AI.\I.,G  VOLTAGE  kaTIo 

«2= 1 0 . * » l«A ISOl/20 , )  JFrAnSFEPEo  vulTAGE  kATIu 

nUH  I  l  I  *t  A  :4  S  A  1  f  Pu*E 
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UO  39 
0040 
004  1 
0  0  **2 

0043 

0044 
0045 
004b 
0047 
00  46 
u049 
0050 
0051 
0052 


0053 
0054 
0055 
005o 
0057 
U05d 

0059 

oOoo 
oOo  1 
00o2 
ooo3 

0004 

o0b5 

OOOO 
Ooo7 
OOoO 
0069 
0070 
0071 
0072 
0073 
00  74 

[ 

I 


1 

[ 

I 


ANGrtlsCMPLXld.,0.) 

A.\.GvT=CMPLX(0.,PhITV7 
£H1Ps£1H*CEXP (ANGhT) 

£VTPs£1v*C£XP(AnGVT) 

C  HOPIZUNTAl  TPAnSaiIT  CUMPONEnI  Air h  ANTET.nA  X-CUOPlING 

EHT=EmIP*P1*EVTP»P2 

C  vEKTICAL  THAmSiAIT  COMPONENT  AllM  AnTENna  X-CUUPlING 

£vi=£vrp*Ki+tHrp«i<2 

PMlPFi=o.  IPmASE  SnlFT  1 0  DECEIVED  MuP  SIG 

PMINV=o.  IPmASE  ShIFT  TO  PtCElVEu  VtPT  SIG 

TYPE  6003 
ACCEPT  6  002 ,  NS  TEP 
TYPE  6001 
ACCEPT  t»U02,nFFT 
11111  CALL  Pump) 

TYPE  6000 , I  PUN 
L 

c  ***  all  x ivpu r  oata  is  peao  fpqm  pill  Input  ro  ifile 

c 

ACCEPT  bol2,lFILt 

CALL  ASSIGN  122.  XFILE.o, 'koo*  J 

pEAO  122,  o0o2  J  No  ISE  IEnTE*  o  F  UP  ..OISE  OFF,  l  FOP  NOISE  uu 

PtAO 122 , 0004 J P IFdA  JPtCtlvEh  IF  0Anu*1O1H  1 n  mEnTZ 
k£ AO (22 ,6004)RNFoB  IPECEIVEP  NOISE  FIGUkE  In  Ob 

PiVF s  1 0  .  «*  (KivFOd/20 .  J  JpECEIVEk  *01  St  Fl&OPt 

C  VAKIANCE  s  nx T 8nF 

VA«Is(50.»1 .38E-23«290.«PIFB6*«NF)  J50  Umk  IPPEOAnCE 

SO=S«XPTlWA«X)  SSTAMJAKO  DEVIATION  s  SOHT  IVARIAmCE  J 

kEAO (22,6oU4ICF  X TRANSMITTER  CENTER  FREQUENCY 

KcAOl22.o004JFd«  1F«EJUEUCY  AGILITY  8AND*IuTrt 

Kt  AU  122 , 6o  o  4  )  PnF  XTkAivSmIT  PulSE  pEP  FpEO 

uF=F6A/FluAT  l NsTcP- 1 J 

ML  AD  ( 22 ,  O0u4  J  P  1  PrtK  }AvEk*GE  XmIT  Pam  Ml.PIZ  OP  YtMT  CHANNEL 
L  NhEN  XPITTEP  TUmVEO  OA  ll.E»  XMIT  Poalk/£) 

mEAU(22,o00  4JCK  iCO«PPESSI on  maTIO 

p6ADi22,o0o4;GaInA  jai»TE”vAA  GaIn  In  Oo 
anTG= 1 0 . **  IGa Ina/ l o . j  1 an TENna  gain 

ANtG2sANTG»*2. 

m£  Au  122, 600  4  J  P  A'lOt  'RANGE  In  .mETEpS  To  CELL  OF  If.  TEkEST 
PAnGE4sPanGE**4, 

n£ AO  122,0004} DBLoSS  JSYSTEm  LOSSES  IN  DO  fu»  lIThEP  m  OP  V  CnANNtL 
5LUSS=10.**(UdLUSS71o.)  I  SYS  TEM  LOSSES  EITmEm  CHANNEL  M  op  V 

kEAD  122,6004) ASCAlE  ! AMPLITUDE  SCAcE  max  SCAlE 

C 

C  Trt£  NEXT  LINE  MUST  at  NjodtP  OF  SCATTEmS  To  oE  OfcAU 

C  F  POM  input  fill 

L 

C  TmEn  t AC m  SuCCttDI'.o  lInE  \ILL  CMApaCIEmIZE  ThE  SC*TT£m£mS  as  FULLUao 

c 

c 

C  SCATEPli, l)=rrPt 

c  cnIE*  i  FOP  Flat  Plate 

c  2  FOm  DImEDPAL 

C  3  FOP  IpImEOmAl 
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0o73 
00  7  0 
007  7 
0070 
0079 
0000 
0061 
0000 
ooo3 

0009 

Ooo5 

0006 

U0o7 

UUO0 

oOo-* 
oo9u 
uo91 
Ou90 
0093 
0U99 
9995 
UU9b 
0097 
O09B 
0  0  99 
u  1  00 
olol 
0  1  0  0 
0 1 0  3 
O  l  0  9 

0io5 
o  loo 
0  l  o  7 
o  loo 
0  1  o9 
olio 
0110 
olli 
Ol  19 
0115 
olio 
oil? 
olio 
Ol  19 
0100 
0  1  0  1 
0100 
0103 


30 


60 


9  F OK  0 IPOLE 

SCaI£k(I,0)=  SI0E  (SO  PETERS) 

SCATEK(I,3)  =  UKItNlATIUN  ANGLE  IN  DEGKEES 

SCA  IEnU,4)=0.vE  *AY  OlSlANCE  FROM  LEADING  Edge  OF  RANGE  CELL  (METERS) 
KEAO  (00,6OO0)NSCAT 

call  plotioj 

CALL  V 1 9CS0 ( 4) 

T  TPt  6010, IF ILE, ASCALE, NSCAT 
DO  5o  1=1, NSCAT 

KEAO  100,6011) (SCATtK(l,K),K=l, 4) 
lYPE  6009, (SCAIEH(I.K) ,K=1,4) 

Call  V19CSZ11) 
call  close  100) 

DO  bo  1  =  1, NSC  A  1 

SCaTER(  I, 3)=SCAltK( l, 3) «CDK 

IFoIF'lsl  (SI  ART  FR£u  a)  STEP  1  OF  UP  RAMP 

ASNRm=u.  (Initiate  ACLUMolaTu*  FOk  NDlSt-FK£E  H  Chan.EL  SIonaL 
aSnk  v  =o . 

DO  000  I F  =  1 ,  N  S  T  E  P 
F  =F Kt H l IF  ) 

LAmOa=C/F 
PH1F=P14/LAMUA 
£maCC=CMPLX(O.,0.) 

CVACC=C9PLX(u.,U.) 

00  lou  1=1, NSCAT 
PhIFU=Ph1F»SCAT£k(I,4) 
call  GETSM(I) 

ANGH=CNPLA(0.,Prtl»h9pHlFO) 

A'.oOsC-PL*  (o. , PH iKV  +  Prt I FO) 
t  rmlKStK  T  «CE  XP  (-a  ,GH ) 
tTvENT=EoT*CcXPl-ANGV) 

SC ALE =K SC ALE (LA*OA  J 

tH=SCALt* (S^aTaX 1 1 , 1 ) *E  T NOR  ASM A  TN  X l l,0)»ETvERf) 

EV  =  SC AlEaISpA UX  (0, 1) *ETkOh»oi*ATk« 10, 0) »tT VtA  T) 
trtlP=EM  (SAVE  PUKE  RECEIVE  hur  S1G.»al 


(INITAITE 

(Initiate 


ACCUMULA  TOR 
ACCUMULATOk 


(EFFECTIVE  HUR  PHASE  SHIFT 
(EFFECTIVE  VERT  PHASE  SHIFT 


tH=EH»nl *£v»R0 
Ev=Ev»Rl*EHTP»K0 
tMACC=£HACC*EH 
10O  tvACC=EVACC*tV 

IF (No IS£«tD«0)u0T0  100 
vn1=VNUIS£ (SD) 

CALL  wANPHlPvlJ 

V  N0  =  V  ID  I  St  (Sj) 

CALL  <»iFn(F(J) 

VNhI =vvl *CDS IP ,1 1 
VNH  i  =  V  «  1  »SI (P  ,1) 
vi»vI=v,0*Coa(P  ,0) 

V  >.  v  i,  =  V  N  0  «  S 1  *  l  P  i  <  0  ) 
VNHSC'PLX  IVCHl  ,  V  -|Mil) 
V«V=CiPL<(y  .  V l »  VnVn ) 

AS NRH= A SnKH+CAoS lEHACC) 
AShKVSaSnH  V ACAoS lc  V ACC ) 


liSE  PdRE  HCK  SIG 


I  AC C u«ll *  T E  H 
(ACCoMDLATt  v 


■  .Ulst-FHtt  SIGNAL 

Nl,lSE-FKEt  SiG«AL 


T6 


FuKTkAn  lv 

V01C-03F*  TnU  *8”UCT”8*  UOSoSSOd 
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01*4 

EhaClsEhaCC+ONh  IAui)  UUISE 

TU 

H  SIGNAL 

01*5 

tyACC=tVACC»VNy  1  toil  \G  I  a£ 

TU 

y  SIGNAL 

Oltb 

1*0 

LUivT  1MU£ 

o  1  *  7 

hIsKEal  IEhaCC ) 

0  1*8 

n«  =  Al  MAGUnACCJ 

01*9 

y IsHEALlEOACC) 

OliO 

VU  =  A  1MAGIEVACC J 

olil 

HZA.Vt»  =  U. 

oli* 

y  T  ANuso  . 

0153 

IFlHl.EJ.O.U.ANU.Hu.t'l.O.OJGulU  130 

0155 

n2A»G4AlAty*(HU*Hl) 

blib 

130 

IFlvI.tJ.O.O.AwD.vD.EU.O.OJGulCi  140 

0138 

vTANG  =  ATAAl2iyuiylJ 

0159 

14o 

dE  1  »  =  VT  AnG-hZa.<G 

oi«u 

o  1  4* 
U14  i 
0145 
01**7 
ol4d 
0149 
0150 
0  151 

015* 

0  153 
0154 
0155 
u  1  3o 
u  1  5  7 
O  1  58 
0  159 

0  1  OO 
0  1  o  1 
Olo* 
0  1  bJ 
0  105 
0167 
0  1  66 
0  1  69 
0  170 
017* 
ol  73 
ul  74 


L 

c  ar  OEMivlUO*:  aEiA  Is  ZEaoiol  it-  either  hohIZOntal 

C  OK  y£*TICAL  ANGLt  IS  ZEKU. 

c 

IF  InZA'OG.Eu.O.O.UK.  V!  ANG , EU . 0 . 0 J 0E T A =0 . 

SET  As AnUO  (tfETA.PI*) 

IF  lbtTA.oT.Pl Jo£TA  =  0£  U-PI* 

IF  iritTA.LT.-PI)a£TA=Kl<;*0£!A 

sat  T  AH  l  IF  )  shZA.<G 

SbtrAVlIFJsVTAAiG 

SF  kEu  C  IF  1  sF 

HliSAvlIFliMU 

*ISAVUF}=HI 

VUS»VUfis*a 

vIsAwUFlsvI 

AHSAy(IF)sCAdSlErlACC)  'CALCULATING  PEA*  tiUHlZ  AMPLITUDE 

AVSAv(IF)=CAoSltVACC)  iC ALCULA I  1*6  PfcAA  V£kT  AMPLITUDE 

at  I Sa  v ( IF ) =  d£  T  A 

kmCSAO  l  IF  )  =  v  105***C''PLA  I  (nl»vol  ,  lH«»y  1 J  ) 
lhC  Sa  y  I  IF  )  =  V  lus**»C  "Pt  A  llnI»VLJ  >  l  n  «  ♦  y  1 )  I 
LOtT*  =  ATAN<i(l«o*»IJ#lttI-Vu>) 

KtttTAsATA-yailHu-yi),  in  I  »vu)  ) 

CBcTA=LdETA-K0tTA 
Lot  T  A  =  A  100 (Cs£ TA,K1*) 

IF  lCotTA.GT.pl  )CbETAsCdtTA-Pi 
lFICoEt  A.LT.-PI)CoETA=taETA«.Kl 
CofcTASUF)=CoETA 
*oo  CO«ri.«UE 

nsaLuGIOCFlUAU.NFFT  I  J/ALOGlUl*.7+0.5 

IF l I kUn.Ej. 1 ) C ALL  ANUKM  Hi,  i«FF  T  ,  NS  TEP ,  ri  IG) 

5NHM=99.49  IVALUE  IF  NU1SE  IS  luKNtD  UFF 
SNk  y  =  9y . 99 

IF  Ou  I  s£  .£0 . 0  )  uQl  0  *10 
C 

C  C  ALCuL  A  I  £  AyEHAoE  9mS  lUlSE-FKEt  SINGLE  PuLSE  SnK  Fu4  EACH  LHA'.nEL 


0176 

0177 

C 

o!7«  *to 


S9KH  =  *o.«AL0Glo(»lU5K**ASNKH/'(FLUAUi<S(tP)»SU)J 
S9K  V:*U  .  *  *Lt>t»  1  O  (  V  1  US6*«  AS  in  y  /  (FLo*t  tiySltt')*5D)  ) 

CALL  SnaIISaI 
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o  lit 
olol 
0  1  0  3 
01o« 

o  1  oS 
0  lob 

ola7 
o  1  Ort 
0  loS 

0  iso 
o  1  S 1 
OlSP 
o  1  si 
OlSO 
U  1  S3 
0  1  So 

0  111 
ulSb 
o  1SS 
uPoo 


OPol 
OPoP 
oPul 
0pu3 
OPu) 
OPoS 
oP  l  0 
oPl  l 
oPl3 
oPl« 
Op  1  3 
OP  l  o 
OPl  7 
oPia 
op  1  s 
oPeu 
Opel 
OPPP 
Opel 

oPps 
oPeS 
oPpo 
oPP  / 


OPeo 

oPps 


if  iiA.sousi»,p)  ,tj.p)icpr=i 

iFUA*OUS.*,l).E„.0)».010  P50 
CALL  AaSlG*  tb.'LP:  ’.0) 

API  It  Ib.OOOO) 
uu  Ppo  isi.ssrtP 

PPO  API  It  lb,bu07)sFPt  Jill  .rllSAOU)  ,  HU$AV(I  )  ,  OlSAVlI  J  ,  VuSAVl I) , 

1  AnSAvll),AVsAYlI),dE rSAOll) ,SdE  r»Mll J r SBETAv  UJ 

call  close  to) 

PSo  call  ASSIGMPP,  'UKSPISSlM.NAM  ',0,'PoO') 

IC1W=0 

300  CALL  PLUflOi 
if  imsm 

X*-Ii*=(CF-Fap/P.)/l.t9 

x.sa*siCF»Foa/P.)/1.ES 

YSiMSO. 

IMA  X  =  ASCAL£ 

Call  AX£S(A,XHiM,XrtAX, IXMAX-XPIn)/50.,YMIN,YMAX, 

1  (»MAX-TMIN)/SO.,XMlN,YMlI»,d40.,b0O.,15o.,10O.) 

pE*oipp,oOoo)sxaxI3 

AtAljl^2,oOOo}i\YAXlS 

CALL  LAritLlA,SXAXl3,.YYAXIS) 

CALL  ntAOER 
C 

C  PLJTTlMG  nORIZ  OR  RnC  AMPLITUDE 
C 

UO  sou  IF=1,rSTEP 
x=SFkEu(IF)/1.E9 

lFUCIR.EiJ.OlYSArtSAVUF)  JGtTUbG  P£AX  HURIZ  AMPLITUDE 
iFlICIP.tQ.l)Y=CABSlRHCSAY(IFJ)  IUETTInG  PtAA  PhC  AMFLITuOt 

IF  l  IF  IkSI  .E'J.IICALL  LI«E(A,*,Y/uVOLTS,o) 
lFIRSr=0 

SOO  CALL  List  l  A,  X,  Y/uy.JLlS,  1) 

IF  UCPT.t’J.O)  oOYO  s  1  0 
CALL  RPOLPY 
CALL  STALL 
GulG  spo 

<4 1  u  ACCEPT  boOb.IA.iS 
<*Po  CALL  PLOIlo) 

IF i»ST=l 

XMliv=lCF-FbA/P.)/l.tS 

XMAX=tCF»FoA/P.J/l.E9 

Y«ISSO. 

YMAX=ASCALE 

Call  axes  1  a,  X.«ln,  xmax,  ixmAx-xmin) /So. , Y'U  ymax, 

1  (YSAX-YMl,l)/50.,XMllV,r*I'*»b‘*U.,b00.»130,,100.) 

kEao IPP.bOob)  SxAX'S 
st al VPP t ooOb )  .YAX13 
call  LAOtL  C  A  ,  ||X  »  X  lo  .  ..Y  A  X  IS) 

CALL  rttAoER 
C 

C  PLOTTl.b  wtKT  dp  L"C  A-plITuoE 
C 

00  SoO  IP  =  I»iiSlEF 
X  =  sF  *fc u 1 1 F ) / 1 . t9 
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0230  IF  l  IC  IK.ED.U J YsAvSAV  l  IF)  1  GET  TING  Pt An  VEkT  AMPLITUDE 

U232  IF  (1CTK.ED.  1)  YsCABsUHCSAVlIFi;  ILUAdInG  PEAr  l«C  AMPLITUDE 

0234  ifufinst.ej.dcall  lin£(a,x,y/uvolts,o) 

o23b  IF  IkS>  TsO 

0 237  buy  call  l!n£  IA, A > T /uVulTS# 1) 

0238  IF(ICPY.EQ.O)  GO  1 0  510 

o24u  call  hROCPY 

o24l  CALL  stall 

0242  GUIU  320 

0243  310  ACCEPT  600b«lAivS 

0244  32o  lFIPSTsl 

u243  XMlN=(CF-FdA/2.)/l.£9 

0240  XMAXs  CCF*FdA/2.  J  /  1  .£<* 

U24  7  TMIN  =  -lau. 

0240  YMAX=1B0. 

0249  call  PLOrtO) 

o2bo  CALL  AXESCA. XMIN, XMAX, (XMAX-XM1N) /50., Y!<ilN, YMAX, 

1  l YMAX- YM IN)  /SO.  .  XMIN,  YMlN,  84  0.  ,bu(i.  ,  150.  ,  l  0  0.  I 
0231  HEAD  122, oOUb) NXAXlS 

02b2  K£A0l22.o00b)NYAXIS 

0233  CALL  LAoELCA, NaAXIS, NYAXIS) 

o234  CALL  hEADEk 

C 

C  PLOTTING  PHASE  b£  TmEEN  MOV  OR  RHC  &  LHC 
C 

0233  DO  600  IFsl.riSTEP 

o23b  X3SFhEW(IF)/I.£9 

0237  IF IICIK.EQ.O) Y=BETSAV  t IF)»180./PI 

0239  IFUClK.EQ.IJY=CoETAS(IF)*ldO./PI 

o2bl  IFCIFIkST.EO.DCaLL  LINE(A,X,Y,0) 

02o  3  1FIPSTS0 

o2o4  bOO  CAlL  LINE  (  A.  X,  Y  ,  l  ) 

02b3  IFlICPY.tl.OJ  GOTO  610 

o2o7  CALL  mkOCPY 

o2oa  call  stall 

o2b9  GOTO  o20 

02/0  6 1 0  ACCEPT  6uOb»IAiiS 

u27  1  o2o  COivT  1'VUE 

0272  DO  1500  lFFTsl,2 

o27  3  DO  lolo  I=1,nFFT 

o2 7 4  1010  VALUtlU=CMPLXlO.,O.J  12ER0  OUT  COMPLEX  dJFFEH  VALUE 

0273  DO  1020  I  =  l,.vSTEP 

02  76  IF  (ICIh.EG.O.A  vu.IFFl  .E0.17  VALUE  IU  = 

1  CMPLX  InISA  v  (  I )  r  HUSA  V  1 1)  1  JLUAu  BUFFER  *ITf>  fiuHIZ  I  4  k 

0  2  78  IF lIClrt.EU.O. AND. IFFT. EU.2) VALUE  II  )  = 

1  C'HPL  x  l  AHSa  V  (  I  )  »  V  .  )  ILuAU  O  offer  REAL  PakT  <«ITH  hURIZ  A  mP 

o200  IF  l  IC  IP. ED.  I  .  A.yD.  IFFT  .ED.  1  )  VALuE  1 1 1  =RhCSa  y  (  I )  ILUAU  BUFFER  .iITn  ■<" 

u2o2  IFlIClK.EJ.l.  A|*0.  IFFT.  ci.2)VALUE  ll)  = 

1  CiAPlXIlA3SIHHCSaVII))»0.)  ILUAU  BUFFER  REAL  PART  aITh  NHC  PEaa 

0204  1020  COnTIvuE 

o2o5  Avgval»CmPlxio.,o.I 

020b  GO  1 0  l 0  7  o 

C 

c  Rf'uvE  uc  value  from  ff  t  input  uf  implemented) 
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C 

C  00  1050  I  =  l*,*srEP 

L 1 050  AVGVAL=AVGVAL»VALUdli) 

C  A VGV al=A VGV AL/FLUA f ( wS IEP ) 

C  00  1060  1  =  1 ,  nSTEP 

Cloou  VALUtliJsVALuEdJ-AVGVAL 
u267  1070  COiv  T  1'VUE 

C  CALL  OHnATE  (VALUE, iv)  1FFT  INPUT  /itIGHTING  IF  C  ALL t  U 

02o8  CALL  iLOGiM IN,  VALUE,  *1.1 

0289  CALL  PlOI  (o) 

0290  0ELX  =  C/(2.*Ft3rtJ 

0291  1XMAX=IF IX (DELX* (NFFT-1 )+0.5) 

0292  INfcMAN=M(J0(IXMAX,5) 

0293  IF UKEMAN.tO.OJGOTU  1060 

0295  IXMAX=IXMAX+(5-IKEMAN) 

0296  1 060  XMIN=0. 

0297  XMAX=FLOAT (IXMAX 1 

o29d  Yn1N=u. 

0299  Yl“AX=l. 

0  30  0  call  AXES  t  A,  X.9lN,  XMAX,  (X  «AX-XPilN) /So.  , 

1  i,ii.4,r9Ax,cr'<Ax-rMiN)/5o.,xMiN,rKiNId9u.,boo.,i5o.,ioo.) 
o3ol  tit  Au  (22*  600b)  NX  AX  I S 

u3o2  nEAO  (22,b006).*YAXlS 

o 3o3  call  la6Elia,nxaxis.ntaxisj 

03u8  CALL  ritAOtK 

C 

C  PLOTTING  FFT  OF  HOKIZ  OK  RrtC  CHANNEL 
C 

0305  00  1100  1  =  1*  "IFF  T 

0  306  X=ll-1 J*UELX 

vSu  7  Y=CAoS(VALUE(I))/PlG 

0  3  0  5  call  L l Vt ( A , X , 0 . , 0  I 

o3i*9  lloO  CALL  L I vt ( A , X  * Y *  1 ) 

0  310  IF  (  ICPY.t  9.0)  GOtlJ  1110 

0312  CALL  "KOCPY 

u  3 1 3  CAlL  SIALL 

u31<*  GO  TO  1120 

0315  1110  ACCEPT  0006, XAnS 

o3ib  ii2o  call  plot  ioj 

0317  XPINSO. 

0  3  Id  XMAX=FlOAT  (  1  XiXA  x  j 

0319  Y9XN=-l«0. 

u  32>>  YHAX=160. 

1*321  CALL  AXES  (A,  XWlN,  XMAX,  (X*  A  X  -  XM 1 1*  )  /  50  .  , 

1  YHlN*Y*«AX,(Y9AX-YMlN)/5O.,XMIN,YPH,rt<40.,bO0.,15O.,lO0.) 
o 322  kEAO  (22,60o6)nXax  is 

o323  KEAO  ( 22  * 6u Oe 1 n Y A x X S 

0329  CALL  LA9ELU,9XAAIS,.*YA<16) 

o325  CALL  hEAOEk 

C 

c  plotting  homu  ok  phc  fft  phase  a.*gle  data 

L 

0326  00  1200  X  =  l,ivFFT 

0327  X=U-U»UtLX 
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u32H 
u3i:<> 
0350 
0332 
o334 
0353 
Ui3b 
os3i 
</330 
0340 
0341 
0342 
0345 
0  344 
o345 
0  34b 
034/ 
0340 

oS30 

0  332 
oS34 

033b 

0337 

0350 


0  350 

0300 
0  3o  t 
03b2 
O3o3 
0364 
0303 

0360 

0  3bT 
O3b0 
o  3o4 
o3/0 


u  37  1 
u  57  2 


YSATA.,2(AIhAGlVALUtUl  >  , RE AL t V ALUE C I )  )  1 
Y=AmuDIY,PI2) 

IFIY.GT.P11YSY-PI2 
If CY.LT.-PI)Y=PI2*T 
y=r*too./Pi 

CALL  HNE(A,X,G.,OJ 
1200  CALL  LINtlA,X,Y,ll 

IHICPY.E0.01  GOTO  1210 
CALL  HROCPY 
CALL  S I  ALL 
GOU)  1300 

1210  ACCEPT  600b, 1 AnS 
1300  C(7ATi»U£ 

00  2300  lFFTsi,2 
00  2010  I* 1 , RFF  T 

2010  VALUEtIJ=CMPLX(0.,  0.1  UfcHQ  OUT  COMPLEX  t»Ufft«  VALUE 

00  2020  I  =  1,.VSTEP 

IFUCIR.EU.O.ANO.IFFT  .E0.1)VALUEU1  = 
l  CMPlXIVISAVIU.VOSAVU))  1LUA0  BUFFER  "ITh  v£RT  XAU 
IFUClR.EQ.O.ANO.IF  f  l  .  1 3.2)  VALUE  Ul  = 

1  Cm?LXIAV3AVII),0,J  ILoAU  BUFFER  REAL  PART  *1TH  VERT  AmP 

lFaCl4.EU.l.A^0.IFFT.fcQ.l)VALUtin=LHCSAYH)  1L0A0  BUFFfcR  *ITh  L* 

IF  UCl*.  EG.  l.A.vO.lFFT.  £3.2)  VALUE  111* 

l  CmPL  X  CC  A8S ILHCSA V  C 1)1,0.)  ILOAO  duFFEK  REAL  PART  blTM  L«C  AMP 

2020  CONTINUE 

AVGVALsCMPLXtO.,0.) 

GOTO  2070 
C 

C  REMOVE  OC  VALUE  FROM  FFT  INPUT  UF  IMPLEM£M£0) 

C 

C  00  2o30  lsl,.v3TEP 
C2o3o  avom al=A vGval»VAlu£( 1 7 
C  aVGV ALsAVGVAL/FLUAT (NSTEP) 

l  JO  2ob0  Isl/.lSTCP 
C206O  VALuEll)=VALUEUJ-AVGVAL 
2070  CONTINUE 

C  CALL  OfiAAfEl  VALUE,. V)  IFF  T  INPUT  aEI&mTI.hG  IF  CALLEU 

CALL  NLOGVtN, VALuE, tl.) 

CALL  PLOT  tOJ 
amIusO. 

AMAXsFLOAT IIXMAX) 

YMlNSO. 

YMAXsl, 

CALL  AXES  (A, XrtlN, XMAt , UPAX-XHIN) 730., 

1  YMt.l,  YMAX.  tYPAX-YMlNj/50  XMIN,  YMN,  840.  ,600.,  130,  ,  100.1 
READ  (22,t>00«)'VXAXi3 
rlEAO  C22«6UUb).xYAX  1$ 

CALL  LA6EL(A,4XAXIS,NYAXISJ 
CALL  MtAutb 
C 

C  PLOTTING  FF1  OF  VERT  JR  LMC  CMNL 
C 

00  2100  1*1,  ,FFT 
XSU-1)*1)>-LX 
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PACE  009 


0373  Y=CAb3lVALLE(.IJi/'uiO 

o374  CALL  LlVtUil.D.iU) 

v>  3  7  3  EtOO  CAul  L  live  (A  ,  A,  f ,  1 J 

u37b  IFUCPy.Eg.OJ  0010  2110 

o37<l  CALL  HKDCPY 

0379  CAlL  STALL 

o56o  C0 1 0  2120 

i>3o  1  2110  ACCEPT  bu0b,IA*5 

uio2  2120  CALL  PLOT  10) 

o3bi  XHIftSO. 

u3o<*  *nax=floatiixmax; 

03<s5  TMl,v  =  “l8v. 

oibb  YMAX  =  lavl 

USb7  call  AXES  (A.XmIN,  XMaX,  (X,«a K-Xhl:<i)/bU., 

I  YMllY,Yi*AX,CYMAX'»Y*ClV)/’S0.,XMl.<,Y*IN,84i/,,bOO,,130,,19O.) 
U3S8  kEAU  l2a.bu0b)NXAXIS 

0ib9  *E  AO  iaa,bO0b)*YAXIS 

oi90  CALL  LASEHA.NXAaIS.imYAXIS) 

u391  CALL  heaoEk 

C 

c  plotting  vEht  o«  lhc  ffi  phase  angle 
c 

oiva  uo  2200  i = l ,  ™f  f  t 

Ui93  XSU-1J«0ELX 

0  394  YSATAVa  IA IMA6C VALUE ll)),«eAL(VALuE<IJ i) 

1>39S  YSAMUiHY,P12) 

039b  IFIY.ST.PI)Y  =  Y-P12 

0398  1F(Y.LT.-P1)Y=PI2AY 

i>400  Y5Y«iau./PI 

04g  l  CALL  HN£tA,X,U.tO) 

04v2  aaoo  CALL  Ll'vE  (A  ,  X  ,  Y  ,  1  ) 
u<H)J  iP'UCPY.tG.O)  bl)  TO  aaiv 

v«w5  Call  hhucpy 

u«ub  call  stall 

0<*U7  OOIO  2300 

U«U«  2210  ACCEPT  b  0  V  b  ,  i  A  H  S 

u“u9  aSuo  continue 

U<nu  IF  UCIP.tO.DGuTO  2bo0 

vsia  ici«=i 

U«»13  GOTO  300 

ua i 4  atoo  call  close  caaj 

0415  ixUtvalKUNtl 

041b  IF  1 iKUn.tQ.ajGoTo  30UU 

uoia  ooTo  niu 

U«19  3  0  UO  CALL  V14CSZ14) 

C 

C 

uvao  boia  FOKMATisAaj 
u«al  oOll  F  OrtrtA  T  1 4F  2u . 0 ) 

«4aa  boio  fokma r ( i x» *f ile  na*e:  '«a2x 

1  lx, 'SCALE!  '  ,r  7 .0/ IX  , '..LPtsck  SCATTExs:  »,le/) 
uua3  bOQ9  FOkhAT C) X,31F7.a, ' , ' ) ,F7 ,ai 
04<4  b008  FOH.AAU  1  X,  T5,  'fkEU'  , 

1  118,'HURZ  1',131,'HOKZ  0',T44,'VEhT  I«,Ts?,'VtHl  O', 
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PAUE  010 


0423 

6007 

0426 

o006 

042/ 

60OS 

0428 

6004 

0429 

o0u3 

0430 

6002 

0431 

6001 

0432 

60o0 

0433 

1  Uo.’HOKZ  AMP* , T83,  'VERT  AMP', 196, 'BETA', 

1  T 1 09 » '  rt  BETa',T122, 'V  BETA') 

FORMAT! lv!2X, 1 PE 11. 4)) 

Format (2ua2) 

FORMAT!'  ANTENNA  ISOLATION  IN  OB') 

FORMAT !F 1 0 . 0) 

FORMAT ( •  NUMBER  OF  FHEUutNCr  RAMP  STEPS?'/) 

FORMAT l 16) 

FORMAT !  / '  NUMBER  OF  FFT  POINTS  IlESS  Ok  EJUAL  256'/) 

formatc  run  number  ',13//'  uata  file  for  scatterers*/) 

ENO 


FOKTkAP. 

IV 

SIUHAGt  MAP 

N|  Atot 

OF  F  St  1 

A  T  Tk ISO  TES 

MAX  IS 

OOOOUfe 

i .4  T  £GEH*P 

AHHAT  IPO) 

N  T  A  X  I  S 

00005b 

Ik  ItGEH*P 

ahka r  IPO) 

A 

UOoleb 

HE AL  *4 

AkkAY  IPu) 

Ldt  T  A 

001200 

HEAL*4 

V  Ak 1  ABLE 

LAMUa 

uoieop 

HEAL  *  4 

vahiasle 

AVbVXL 

VO  12  lu 

COMPLEX  *  0 

VAklASLt 

tH 

OOlPdO 

COMPLEX*** 

VAHIASlE 

tv 

OOldiO 

COmPlE  X  *0 

VAHIaSlE 

t  frlUH 

0O1P4O 

COMPLEX*a 

vahiasle 

ETVEkT 

ooipso 

COmPlEX«S 

VAHIASLE 

AivIIjM 

V  0  1  P  o  0 

C  OrtPLfc X  *  S 

VAHIASLE 

AkbV 

voia/v 

COMPlEXaS 

VAHIASLE 

tMTP 

V01 JVO 

C UmPl£X*S 

VARIABLE 

tv  TP 

001310 

COmPlE  x  *a 

vahiasle 

tHT 

VO l 3PO 

COMPlE X*d 

V  AN  I  AdLE 

tv  I 

ooiiio 

COMPLEX  *tt 

vahiasle 

AkGH  T 

VOI 340 

C  umPlE  X  *  S 

vahiasle 

Ai\bV  1 

001 55o 

COmPlE  X  «  d 

V  Ak  I  ASlE 

CMACC 

oO 1  iso 

C  UmPlE  X  *o 

VAN  IxdLE 

L  V  ACC 

oOli/o 

complex *d 

VAHIASLE 

VNn 

OU 14u0 

COmPlE  x  *  s 

vahiasle 

VMV 

001410 

CumPlEx*s 

vahiasle 

VIUSH2 

O014PO 

HE  AL*4 

vahiasle 

suwr 

000000 

k  E  Al  *  4 

PROCEOuHE 

UVULTS 

u014£o 

kEAL*4 

variable 

C 

001430 

HEAL*4 

vahiasle 

ICPY 

00  1434 

IHTEGEH*P 

variable 

IkON 

001430 

Ik  ft GEH  *P 

variable 

PI 

001440 

HgAL  «  4 

vahiasle 

i-li 

001444 

k£  AL  *4 

vAhIaolE 

PI4 

O014S0 

At  AL  *  4 

V  Ak I ASlE 

CDk 

00 1 454 

H  t  A  L  *  4 

v  Ak I AslE 

plut 

ooOOOO 

Ht AL*« 

PRuCEUUkE 

V14CS£ 

oOoOoO 

kEAL*4 

PkOCEJokE 

XMlT 

O0000O 

HEAL*« 

PHGCtOuHE 

tIH 

0  0  1  4  S  o 

HEAL  *  » 

VAHIASLE 

tlV 

O0 1 464 

He AL*4 

VAHIASLE 

Phi  TV 

001470 

«E  AL  *  4 

VAHIASLE 

k  1 

O01474 

HE  AL  *  4 

vahiasle 

001500 

xEAL*4 

VAHIASLE 

C«PLX 

OOOOOO 

COMPLtX*d 

PRuCEOuHE 

CtxP 

oOOOOO 

CUmPLEX*S 

PRUCtOuRE 

PrtlHH 

00 1 504 

HE  AL  *  4 

V  Ah  I AdLE 

pH  I  k  V 

ool5lo 

Ht  AL  *  4 

variable 

..FF  T 

001514 

I V T£oEh»P 

V  AH  I AdL  E 

ASSIb* 

oooOoo 

HE  AL  *  4 

PRuCtUUWt 

k  IF  do 

00 1 5 1 o 

KtA(_*4 

v*hI AdLE 

HkF  US 

001 5PP 

HE  AL • 4 

vahiaoLE 

k.«F 

o  0 1 5eo 

k t  AL  *4 

VARIABLE 

V  A"  I 

OolSiP 

kt Au*4 

VAhIASuE 

F 

oo 153b 

HE  AL*4 

V  AH  I AdLE 

FLuaT 

OOOOOO 

he  al  *  4 

PRUCEOUHE 

FORTRAN  IV  storage  map 


NAME 

OFFSET 

AT  TRIBUTES 

ANTG 

001548 

R  EAL  »  4 

VARIABLE 

ASCAlE 

VO  1 546 

REAL*4 

variable 

I 

w  1 5 

IN  t  EGER  »8 

variable 

K 

00l!»54 

INTEGER**! 

variable 

CLOSE 

vOOOOU 

REAL  *4 

PROCEDURE 

ASNNH 

O015S6 

RE  AL  *4 

VARIABLE 

asnrv 

U0l5b8 

n£Al«4 

VARIABLE 

IF 

VO  1 56b 

IN  TEGEN*8 

VARIABLE 

F 

v  0 1 5  7  U 

RtAL*4 

VARlAbLE 

FREQ 

VOOOOO 

REAL*4 

PROCEDURE 

PNIF 

O01574 

KEAL«4 

VARIABLE 

PMIFO 

VO  1600 

REAL*4 

VARIABLE 

getsm 

VOOOOO 

REAL*4 

PRUCEDURE 

scale 

001604 

«EAL»4 

variable 

rscalE 

VOOOOO 

REAL*4 

PROCEDURE 

VN1 

VO  16 1 0 

REAL*4 

variable 

VNOISfc 

VOOOOO 

REAL  *4 

PRUCEDURE 

K  ANP 

_vooooo 

RtAL*4 

PROCEUuRt 

PN 1 

001614 

REAl*4 

VARIABLE 

VN8 

U01680 

RE AL  *  4 

VARIABLE 

PN8 

001684 

REAL»4 

variable 

VNMl 

001630 

REAL«4 

variable 

COS 

VOOVOO 

RE AL»4 

PROCEDURE 

VNHU 

001634 

REAL«4 

VARIABLE 

SIN 

VOOOOO 

«£AL*4 

PROCEDURE 

VNVl 

00164V 

n£AL*4 

VARIABLE 

VNVQ 

V01B44 

R£AL*4 

VARIABLE 

CAOS 

000000 

REAL  *  4 

PROCEDURE 

Ml 

V016S0 

NEAL A4 

variable 

kEaL 

VOVOOV 

«EAL*4 

PRoC£uuR£ 

nO 

0016B4 

n£ Al  *  4 

V  An  I  ABLE 

A  I  MAG 

VOOOOO 

n£AL*4 

PROCEDURE 

Vi 

001660 

REAL*4 

V  AR I  ABLE 

vg 

001664 

HE  AL  *  4 

variable 

MZANG 

00 16  TO 

REALA4 

variable 

y  T  ANG 

00 16  74 

R€Al*« 

variable 

A  T  AN8 

VOOOOO 

rEAL»4 

PROCEDURE 

BETA 

vO  1 7  00 

REAL A4 

VARIABLE 

AMOO 

OOvOvO 

RE AL  *4 

PROCEDURE 

hBET  A 

001704 

r£  AL  *4 

variable 

CBt  T  A 

001710 

REAL*4 

VARIABLE 

N 

001714 

I N  T  £  t»E  R  *  8 

variable 

AL0G10 

VOVOoO 

RE  AL  *  4 

PROCEDURE 

anuRm 

vO VO vO 

HEAL  *4 

PRuCtOURt 

SaR 

VOOOOO 

R£AL*4 

PRuCEUORc 

iS« 

VO  t  7  16 

1n1EGEk*8 

V  An  I  ABLE 

IAN0 

vvOOOU 

1  '<  TEG£R*8 

PRuCEOURE 

ICIR 

001780 

1  yT£titR»8 

VAR1 ABLE 

IF  InST 

vO 1 788 

InT£GEn#8 

van IaBlE 

X  m  I  .y 

V01 784 

RcAL»4 

variable 

XMAX 

001730 

RtAL*4 

V*rI Able 

tmin 

VOl 7 A4 

REALA4 

variable 

85 


FORTRAN 

IV 

STORAGE  MAP 

r.AME 

OFFSET 

AT  TRIdUTtS 

UMAX 

uO l 7  40 

RE AL  *  4 

VARl ABLE 

AXES 

u  0  V  0  0  0 

Rt AL*4 

PRUCEOURE 

LAdtL 

uOuGou 

i.\iteger*2 

PRuCtOuRE 

nt ADtR 

uOoOOO 

RE  AL  *  4 

PROCtOURE 

X 

001744 

KEAL*4 

VARIABLE 

r 

001750 

REAL»4 

V  Ak  I  ABLE 

LInE 

UOoOOO 

IN 1 EGER*2 

PROCtOURE 

MRuCHY 

uOUUuO 

k£AL»4 

PRUCEOURE 

stall 

uOuOoo 

REAL*4 

PRUCtOuRt 

IA.vS 

uO 1 754 

1nIEi>ER»2 

VARIABLE 

IFF  T 

VO  1 756 

InTEG£K*2 

VARIABLE 

\LGGn 

UOOOOO 

In[£GEK*2 

PRUCEOURE 

u£lx 

00 l 7bO 

Rt AL*4 

V  Ak I  able 

IXMAX 

00 1 7b4 

1nTEGER«2 

VAR lABLt 

I F  I  X 

oOvOoO 

InTEGER»2 

PRUCEOURE 

Ht^AM 

o  0  1  7  ob 

IvTEG£R*2 

VARIABLE 

SUu 

uOoOuO 

INTEG£R»<? 

PRuCtOuRE 

COmMUN  BLOCK  /aUkkF/  lEnGTh  uoOu24 


IFSTFQ 

UOOOUO 

1 N  TEGER  *2 

VARIABLE 

luP 

U00002 

1NT£G£R*2 

VARIABLE 

LSfEP 

0  0  0  0  V  4 

1ivTEGER*2 

VARIABLE 

NSItP 

UOOOOb 

I.VTEGER*2 

variable 

OF 

oOoOlO 

RE AL*4 

variable 

CF 

OOuG 1 4 

REAL  *4 

variable 

F  BR 

000020 

kEAL«4 

VARIABLE 

C'*M'*U"l 

-LOCK  /rt*sCT/  uEivG  ( n  ouilJO 

SC«TtR  UI'UUUO  k£AL«4  ARKAr  llvd,4)  vtCTUREC 

S'-’ATrx  uni 1 UU  CUiPLEXab  ARkAT  12,2)  VECTURtl) 

COMMON  BLOCK  /nEAO/  LENGTH  UoOll4 


A  UUL 

oooouo 

i<EAL*4 

VARlAdLE 

-iSCAl 

o0u0u4 

1 iV I EGEk  *2 

VARIABLE 

If  A  in#  A 

VOuOob 

RE AL  »4 

VARIABLE 

ftd  I  St 

UOOO  12 

l,v  TEG£R*2 

VAR  I  ABLE 

k  Andft 

OOuO  1  4 

kEAL«4 

variable 

UtJLLj^S 

U0o020 

k£aL»4 

V  An  1  ABLE 

T 

000024 

I.vT£G£r»2 

AWkAV  12) 

iFiLt 

uOoOiO 

In  I EGEk  *2 

ARR  A  1  IS) 

Su 

«0oo5V 

R£AL»4 

VAR  1  A  B  L  t 

O  I  v» 

O0oo54 

KtAL*4 

VAR  I  ABLE 

5'»N  h 

OOoOoo 

kE  aL  *4 

VARIABLE 

3hKV 

W0o0b4 

R t  Al  *  4 

VARIABLE 

8%hM  l 

OU0070 

K  £  At  *  4 

variable 

So*  r«  j 

o&uo  74 

Rt AL*4 

variable 

S  Alt*  V  1 

OOU 100 

k£AL*4 

variable 

$  vnVu 

000104 

RtAL*4 

VARIABLE 

86 


FORTRAN  IV 


STURAGE  MAP 


NAME  OFFSET  ATTRIBUTES 

SNR  oooiio  real*  a  VARIABLE 


CUMMON 

BLOCK  /ivORK/ 

LENGTH 

0 A2U0U 

MUSA  V 

UOUOOO 

nEAL«A 

array 

125b) 

nISAV 

U020U0 

kEal*a 

ARRAY 

(25b) 

VGSAV 

OOAOOb 

REAL* A 

ARRAY 

(25b) 

V  1 5  A  V 

UOOOOO 

RE AL* A 

ARRAY 

(25b) 

SFREu 

ulOOOO 

REAL*A 

ARRAY 

(25b) 

SHtTAM 

012000 

RE AL* A 

ARRAY 

(25b) 

SbETAV 

U 1 AOOO 

REAL*A 

ARRAY 

(25b) 

oETSAV 

016000 

RtAL*A 

ARRAY 

(25b) 

RHCSAV 

U20000 

COMPLEX  *  B 

ARMAY 

(25b) 

LHCSAV 

U2A0UU 

CUMPLEX*B 

ARRAY 

(25b) 

CBtT  AS 

U30UU0 

R£aL»A 

ARRAY 

(25b) 

AHSAV 

032000 

RE AL* A 

ARRAY 

(25b) 

AVSAV 

O3A0UG 

RE AL* A 

ARRAY 

(25b) 

VALUE 

U300UU 

COMPLEX*B 

ARRAY 

(25b) 

COMMON 

BLOCK  /SIGNAL/ 

LENGTH 

0O0U30 

PTPwH 

000000 

RE AL* A 

VARIABLE 

RAnGEA 

UOOOO A 

REAL* A 

variable 

CR 

000010 

REAL*A 

variable 

ANTG2 

00001 A 

REAL*  A 

variable 

SLOSS 

000020 

RCAL*A 

variable 

P  I  AC 

O0002A 

real*a 

variable 

FORTRAN  IV 
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PAGE  001 


UOul 


0002 

ouui 

O0u4 

0003 
o<>0b 
U0o7 
o  o  u  a 

000  V 
0010 
0011 
00  13 
00  lu 
0013 
00  1 6 
Ool  / 
0018 
0019 
0020 
0021 
0022 
oo  23 
0024 
0023 
0  02b 
0027 
0028 
0024 
0030 
0031 
oO  32 
oo33 
o0  34 
0033 
003b 
o037 
00  38 
oo  39 

oOuU 

oOol 

0042 

0"43 

0044 


SUbRUUTINE  XMII  Ulrt,£lV#Phnv,(VAME) 

C 

C  This  SUOROUTInE  DETERmInES  TH£  I  w  A  NSM  1 1  1  £  0  SIGNAL  POLAkIZaTIUN 
C 

INTEGER  NAME  (27 

INTEGER  KMC (2) ,  LHC 12) , MUR (2) ,VtRl2),nVC2) 

DATA  RriC.LMC, huR, VEK, hv/ ' Rh* , 'C  '.'Ln'.'C  ', 

1  'V£* ,  'R  ' , 'm-' , 'V  '/ 

PI=3. 14159 
l  Call  plOT(O) 

TYPE  10 

10  FORMAT ( '  1  *  RMC 1  / '  2  -  LMC'/'  3  -  HORIZONTAL'/'  4  -  VERTICAL'/ 
1  '  3  -  HORIZONTAL  &  VERTICAL'/) 

ACCEPT  15, IX Ml T 
13  F  Ohm A  T I lb ) 

IFUXMIT.LT.1.0W.IXMIT.GT.5)(,0?0  1 
GOTO  (100,200, 300, 400, 5oOl IXHlT 
luo  tlM=l. 

Elv=l  . 

PHITVs-Pl/2 
NAME ( 1 j =RhC  1 1  ) 

NAmE (2) 2RmC (2) 
kE  T URN 
200  E1m=1. 

E  1  Vs! . 

PMIIV=4PI/2 

nAMEU)sLmC(I) 

NAME(2)=LmC(2) 
return 
300  tlM=l. 

E1V=0. 

PmITv=0. 

NAME  1 1 J  =  Hl)«  ( 1 ) 
iv Am£  (2 )  smOR  ( 2 ) 
kE  T  URN 
400  tlH=0, 
tlVSl. 

PHITVSO. 

nAmE(I)svERU) 

NAME(2)=VER(2) 

RETURN 
500  Elrtsi. 

tlV=l. 

PHITVsO. 
nAME ( 1 ) shV ( 1 J 

name  (2)  snv  (2) 

*E  t  URN 
tNO 


f  ORTRAiM  IV  STORAGE  MAP 


IV  AmE 

OFFSET 

ATTRIBUTES 

iNAnE 

000022 

I N  T EGER  *2 

parameter 

ARRAY  12 

RHC 

O0O024 

1NTEG£k*2 

array  (2) 

LHC 

U00030 

In  TEb£R»2 

ARRAY  (2) 

HOR 

000034 

INTEGER*2 

ARRAY  (2) 

vEr 

000040 

livT£G£M*2 

array  12) 

HV 

O00044 

InTEG£R*2 

ARRAY  12) 

tin 

000014 

KEAL»4 

PARAMETER 

VARIABLE 

ElV 

u0o016 

R£Al*4 

PAKAmE  TER 

VAR  I  ABLE 

PMITV 

oouoao 

m£AL*4 

parameter 

variable 

PI 

000214 

h£AL»4 

VARIABLE 

PLOT 

000000 

k£AL*4 

PROCEDURE 

IXMT 

000220 

InTEGER*2 

VARIABLE 

FORTpAn  IV 
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uou  1 


uooe 

oooi 

oooo 

OOOS 

uoofe 

uuu7 

0008 


SUBROUTINE  6H<vATt (A,nJ 
C 

c  eft  Input  weighting 
c 

COMPLEX  A  In) 
u a  r  a  PiE/b.Ea3ias/ 

UU  1O0  1  =  1  ,N 

.<ATE  =  0.4E3E3-0.4')7S5*CUS(PIE/N»  1 1-1 J)  *0 . 0  79EE*CUS  (  P  It/ i,«  E  «  (  I  -  1  )  ) 
loO  AU)=A(I)*WATE 
p£  T UPN 
LNU 


FORTRAN  IV  STORAGE  MAP 


NAME 

uF  F  St  T 

a  r  tk isutes 

A 

oooo  l  <* 

COiMPlE  X  *  8 

PAPAmE  TEP 

AKRAY  l«) 

N 

00 uu  1  b 

In  TEGEP  *E 

PAp AMt  T  t  n 

V  A  k  I  A  y  L  c. 

PIE 

OuOuEO 

ft£AL«4 

VAPlAdLE 

I 

o0o040 

1NTEGE«*E 

VAPlABLt 

rtA  It 

oouoaa 

PEAC*4 

VARIABLE 

COS 

OOoOoO 

PEAL*4 

PHUCEOUPE 

FORTRAN  IV 


VU1C-03F*  ThU  E8-OCT-8E  oo:o7:0E 
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UUG 1  FUNCTIUN  VNOISE(SO) 

c 

C  THIS  FUNCTION  GENERATES  GauSSIAN  OlsTPIduTEO  NUlSt  vOlTaGE 
C 


UOuE 

3UM  =  U. 

ooo3 

00  lu  1  =  1,  IE 

GO  0*4  lo 

oum  =  oUM-*PAihF  (ijj 

oOOb 

vNulS£=(SuM-b.J  *S0 

ouOb 

return 

oOo  7 

tNl) 

FONT* A* 

IV 

STORAGE  MAP 

NAME 

OFFSET 

ATTRIBUTES 

VMUISE 

0000  lb 

PEAL*4 

VARIABLE 

bL> 

OOoO 1 4 

p£AL*4 

PARAMETER  VARIABLE 

bUM 

oOoOEE 

PEAL«4 

v AP I AbLt 

1 

0  0  OOEb 

l n  T  6GEP*E 

VAp I A8CE 

rt  A  NF 

U  0  u  0  0  u 

PE  AL • 4 

PPOCEOoPE 

U 

uuou3o 

Kt»L*« 

V  An  I  AncE 

AD-A129  502  ANALYTICAL  RESEARCH  BY  COMPUTER  SIMULATION  OF 
I  DEVELOPMENTAL  POLAR1METR IC/ . . (U)  ARMY  MISSILE  COMMAND 

REDSTONE  ARSENAL  AL  ADVANCED  SENSORS  DIR.. 
UNCLASSIFIED  R  F  RUSSELL  ET  AL.  DEC  82  DRSMI -RE-83-7-TR  F/G  17/9 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 
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c 

0001  SUbRUUT InE  RANPH(PHASE) 

c 

C  THIS  SUBROUTINE  GENERATES  UNIFOkMALLT  DISTRIBUTED  PHASE  NOISE 
C 

0002  DATA  PI  /3« 14139/ 

00o3  PHASE-(RANF (U)“.5)*PI 

OO 04  RETURN 

OOub  END 


FORTRAN 

IV 

STORAGE 

MAP 

NAME 

OFFSET 

attributes 

PHASE 

O00014 

RE AL*4 

PARAMETER  variable 

pi 

000016 

REAL** 

variable 

RANF 

OOOOOO 

REAL** 

PROCEDURE 

u 

o0o022 

REAL** 

VARlAdLE 

FORTRAN  IV  V01C-03F*  THU  2B-UCT-*2  OOUtTUQ  PAGE  001 

C 

0001 

C 
C 
C 

U002 

uuoi 

UUU4 

o0u5 


FORTRAN  IV  STORAGE  HAP 

rare  OFFSET  ATTRIBUTES 

RAnF  OOOOtfa  REAL**  VARIABLE 

u  0000 1 4  REAL* A  PARAMETER  VARIABLE 

I  000022  InTEG£R*2  VARlAdLE 

j  000024  1*TEGER*2  VARIABLE 

HAN  oOOOOO  REAL**  PROCEDURE 


Function  ranF(u) 
uniform  number  generator 

DATA  I,J/0,0/ 
RANFsRANU,  J7 

return 

END 
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0001 

C 

L 

c 

c 


0002 

OOoi 

0004 

ooOB 

oOub 

OOU  7 

0008 

UOOV 

0010 

C 

001  1 

0012 

00  t  3 

3 

OOt« 

0015 

5 

oo  lb 

001/ 

oois 

OOIR 

13 

0020 

uodl 

15 

0022 

0023 

on  2  4 

0025 

10 

002b 

uu«:7 

20 

0028 

oois 

o030 

otlil 

23 

oO  32 

0033 

25 

0034 

00  55 

003b 

0037 

SUdkOOl In£  ANORM(N,NFF  T.NSTEP.BIG) 

THIS  SUBROUTINE  DETERMINES  the  BIGGEST  FFT  UUTPUT  F UK  NUKMALlIATIUN 
OF  FFT  PLOTS 

COMPLEX  V ALUE 1 25b }«SMATRX(2»2) 

DIMENSION  SF RED (£56) 

0 1 MENS 1  OH  SBETAHl25b) , SBETAV(25b) 

DIMENSION  SCAT£HU00,4)  ,  AHSAV  125b ) ,  A  VS  A  V  (25*  1  .  BE  T  S  A  V  (25b  ) 
DIMENSION  M0SAVl2So)  ,HISAV(25b)  , VUSAVC2Sb)  ,  VISAVl25bl 
DIMENSION  CBET  AS  125b) 

COMPLEX  KMCSA V ( 25b) ,LHCSAV(25b) 

COMMON  /kOKK/HuSAV, HISAV, VU3A V , V l SA V , SF RE U, SSE T AH, SBE T A V , BE Tb A  V , 
1  RHCSAV,LmCSAV,CBE  T AS, AhSAV, AvSAV, 

l  VALUE 

COMMON  /*KSCT/SCAT£R,SmATRX 

BlG=0. 

DO  3  I=l,NFFT 

value  uj=cmplx(o.,o.j 

00  5  I  s  |  , NS T £P 

VALUE  II  JsCMPLX  tCABSIKHCSAV  1 1)  J  ,  0.0)  IREAL  PAM  aITh  KtiC  AMP 
CALL  VLOGN  lN,vALUt,M.) 

CALL  BIGESTtNFFr.eifiJ 
00  13  lal.NFFT 

V  ALUE II) =CmPLX  (0  . , 0 .  } 

DO  15  I3UNSTEP 

V  aloe C l ) =CMPLX  ICABStLHCSAV ( I J  J ,0,0)  IREAL  PART  «ITH  LHC  AMP 
CALL  NLOGN  (N, VALUE, ♦!.) 

CALL  dlGEST (nFFT,BIS) 

00  Id  I=l,NFFT 
VALOt(l)=C«PLXiO.,0.) 

DU  20  1st, N STEP 

VALUtU IsCMPLX lAnSAV 1 II , 0.1  JKEAl  PAnT  aITH  H0k 1 l  A*p 
CALL  NLOGNIN, VALUE, *1.) 

CALL  dIGEST  CNF FT, BIS) 

DO  23  I=1,nFFT 
VALUE II)=CmPlX(0.,O.) 

DO  25  I*1,NSTEP 

VALUEU1=CMPLXIAVSaVU)  ,0.1  IREAL  PART  rITh  VERT  AMP 

CALL  nlugnin,value,*i.i 

call  B IGEST (nFF  T , B IGl 

return 

END 
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FORTRAN  IV 


stonage  map 


NAME 

OFFSET 

ATTRIBUTES 

N 

00001* 

INTEbEN*2 

pahameten 

VARIABLE 

•\FF  T 

000016 

IvltGEH*2 

PARAMETER 

VAHIABlE 

N5IEP 

000020 

1NTEDEN»2 

PARAMETER 

variable 

bib 

000022 

n£aL** 

PARAMETER 

VAHIABLE 

i 

oouo5o 

lnt  1  EGER*2 

variable 

C**PLX 

OOoOOO 

COMPLEX*# 

PROCEDURE 

CAbS 

OOOOOO 

NEAL** 

pruceoure 

NLuDm 

OOoOOO 

InTEGEN*2 

PRuCEUuNE 

BlbEST 

uOOOoO 

heal** 

PROCEDURE 

COMMON  BLOCK  /NONA/  LENGTH  0*2000 

hOSAv  OOOOOO  NEAL**  ARRAY  (256) 

nlSAV  002000  NEAL**  ARNAY  (256) 

VOSAV  00*000  NEAL**  ARRAY  (256) 

VlSAV  OOoOOO  NEAL**  ARNAY  (256) 

SFnEu  oluOoO  NEAL**  ARNAY  (256) 

Sat  T  AH  012000  NEAL**  ARNAY  (256) 

SBETAV  01*0u0  hEaL**  ARNAY  (256) 

oETSAV  016000  NEAL**  ARnAY  (256) 

khCSAV  o20 000  COMPLEX*®  ARNAY  (256) 

LhCSAV  02*000  COMPLEX*#  ARRAY  (256) 

CdtTAS  O30000  NEAL**  ANNAY  (256) 

AHSAV  0 32000  NEAL**  ARNAY  (256) 

AVSAV  03*000  NEAL**  ARRAY  (256) 

VALUE  036000  COMPLEX**  ANNAY  (256) 

COMMON  BLOCK  / mK5C T /  LENGTH  003140 

SCaTeN  OOOOOO  NEAL**  ARRAY  (loo,*)  VECTORED 
j^ATHX  OU3100  COMPLEX*#  ARRAY  (2,2)  VECTORED 
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PA6E  001 


C 

C 

0001 

00o2 

00o3 

0004 

ooo5 

oOOb 

0007 

uOuB 

0009 


ooio 

c 

0011 

0012 

0013  10 

0015 
uo  lb 


Subroutine  BIGtSI  (m,BI&) 

COMPl EX  VALu£(2Sb),SMATRx(2,2) 

oIHEnSIOn  SFREul25b) 

uImEnSION  SdE  T AH (25b) , SBE T A V (25b ) 

DIMENSION  SCATER(100,4),AH3AV(25b).AvSAV(25b),BETSAV(25b) 

0 ImEnSIUn  HOS A V (25b) .HISAV (25b) . vUSAV (25b) , VISA V (25b) 

DIMENSION  C6E T AS (256 ) 

COhPlEX  RHCSAV(25b),LHC8AV(2Sb) 

COMMON  /«0RK/HUSAV,h1SAV,VQSAV,VI5AV,SFHEU,SBET AH, SBETA V  *  BE TS* V , 
1  RHCSA  V , LHC  SA  V »CBtTAS,AHSAV»AVSAV, 

1  VALUE 

COMMON  / nksc t/sc ater » smatrx 


UO  10  I  =  1,N 
AVALsCAdS(VALUE (I)) 
IF(*vAL.GT.BIG)8IG=AVAL 

return 

ENO 


FORTRAN 

IV 

STORAGE  MAP 

NAME 

offset 

ATTRIBUTES 

N 

000014 

INTEGER*2 

PARAMETER 

variable 

BIG 

OOuOlb 

Rf AL*4 

PARAMETER 

variable 

1 

OQ0034 

1NTEGER*2 

VARIABLE 

AVAL 

oQoO Jb 

RE AL*4 

VARIABLE 

CABS 

uOouuO 

rEA(_»4 

PROCEDURE 

Common 

BLOC*  /aura/ 

LENGTH 

042000 

HQSAV 

OOO0O0 

REAL*4 

ARRAY 

(25b) 

hISAV 

0020o0 

REAL*4 

ARRAY 

(25b) 

VQSAV 

o04ooO 

REALA4 

array 

(25b) 

v  Isa  v 

oObooO 

kEAL*4 

ARRAY 

(25b) 

SFREU 

OlOooO 

REAL  *4 

ARRAY 

(25b) 

SBETaH 

Ol20uO 

hEAL*4 

ARRAY 

(25b) 

SbeTAV 

OlMOoO 

R£AL*4 

ARRAY 

(25b) 

Be  I SA  v 

OIbOoO 

rEAL*4 

ARRAY 

(25b) 

rhCSav 

02o0o0 

CURPLEX*8 

ARRAY 

(25b) 

LHCSAV 

024000 

C OmPL£X*B 

ARRAY 

(25b) 

CBeTaS 

o3o000 

hE  AL • 4 

ARRAY 

(25b) 

AnS*V 

o32oo0 

REAL** 

ARRAY 

(25b) 

AVSAV 

034000 

REAL'S 

ARRAY 

(25b) 

VALUE 

V3o0o0 

COMPLEX*# 

ARRAY 

(25b) 

COm*«uv 

BLOCK  /«** 

ISCT/ 

length 

0o31 40 

scater 

OOoOoO 

rEAL*4 

ARRAY 

(  1  0  0 , 4 

)  VECTuhEO 

SR A TkX 

oo31 oO 

COMPLEX 

•» 

ARRAY 

(2,2) 

VECTORED 

9* 


FORTRAN  IV  V01C-03FN  THU  28-OCT-62  00!07:33 


PACE  001 


OOol 

C 

c 

c 

0002 

U0u3 

0004 

OOwS 

c 

c 

c 

c 

000b 

00o7 

oooe 


function  rscale  (lamoa) 
amplituue  scale  function 

HEAL  LAMOA, IMPED 

COMMON  /SIGNAL/PTPrR,RAN6E4,CK, ANTG2, SL0SS.PI4C 
0 A 1 A  IMPED/50./ 

PRsPTPmR*ANT62*CLAN0A**2.)*CR/ IPI4C*RANG£4*SLUSS) 

PEAK  OUTPUT  VOLTAGE  13  RELATED  TO  AVERAGE  TRANSMITTER 

PORtR  OUTPUT  (mhEn  TRANSMITTER  IS  SNITCHED  ON)  B  T  SURT (2. ) 

RSCALEsSQRT (PR«lMP£0*2.) 

RETURN 

End 


Fortran 

IV 

STORAGE 

MAP 

NAME 

OFFSET 

ATTRIBUTES 

RSCALE 

000016 

REAL«4 

VARIABLE 

LAMOA 

vOOOl 4 

REAL  *4 

PARAMETER  VARIABLE 

IMPED 

O0O022 

R£AL»4 

VARIABLE 

PR 

000026 

real*4* 

variable 

SORT 

OOOOoO 

NEAL«4 

PROCEOUKE 

COMMON 

SLOCK  /SIGNAL/ 

LENGTH  00003U 

PTPRti 

OOOOOO 

R£aL»4 

variable 

RAnGEQ 

000004 

Rt AL»4 

variable 

C« 

OOoOlO 

REAL*4 

variable 

ANTG2 

000014 

REAU*4 

variable 

SLOSS 

000020 

REAL*4 

variable 

PI4C 

000024 

REAL*4 

VARIABLE 
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FORTRAN  IV 


vqic-ojf*  Thu  ae-uCT-aa  oo:d?s3a 


pace  uoi 


OOdl  FUNCTION  FKEUU) 

c 

C  GENERATE  FRtUUEivCT  OUTPUT  STEP  AS  A  FUNCTION  OF 
C  THE  LAST  FHEQUtNC T  RAMP  STEP  TRANSMITTED 

C 


ouoa 

r 

COMMON  /aURKF/IFSTFO, IuP, lstep, nstep, OF, CF,FbR 

Uuu  j 

V 

IF UFSTFQ.NE.UGOTO  S 

OOdS 

IFSTFOSO 

dodb 

FREQ=CF-FBR/d. 

OUUT 

IUP=1 

wove 

LSTEPsl 

uOOV 

RETURN 

uoio 

S 

IF(IUP.NE.l)SOTO  IdO 

ooia 

IF  IlSTEP.NE.nSTEPIGOTO  10 

Oul« 

IUP=d 

uOlS 

RETURN 

Od  lb 

iu 

LST£P=l3TEP*1 

UU 1 7 

FkE J=FRE J*UF 

00  1  B 

RETURN 

UOI  V 

too 

ifilsiep.ne.dgoto  no 

uoai 

IUPS1 

uuBa 

r£  T urn 

Odd  3 

lid 

lstep=lstep-i 

dOdA 

FrEO=FREu-OF 

OUdS 

RETURN 

dddb 

Emu 

FuRTrAN  IV  STORAGE  map 

NAME  OFFSET  ATTRIBUTES 

FREO  uOOdlb  REAL  *  A  VARIABLE 

I  uOudlA  INT£G£R*d  PARAMETER  VARIABLE 


COMMON  BLOCK  /hORRF/  LENGTH  OUOddA 

IFSTFQ  UddOOU  INTEGER**  VARIABLE 
IUP  uOUOud  INTEGER**  VARIABLE 
LSTEP  OOUOOA  INTEGER**  VARIABLE 
nSTEp  uOvdOb  InTEGEr**  VARIABLE 
OF  dOuOlO  REAL  *4  VARIABLE 
CF  OOddlA  REAL*"  VARIABLE 
FBn  OOdO*d  R£Al*A  VARIABLE 
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F OR  TRAN  IV 


V01C-03FV  THU  2B-0CT-S2  00107S4S 


page  ool 


0001  subroutine  GETSM(I) 

c 

c  determine  scatterer  TYPE  and  call  lt*s  matrix 
c 


0002 

COMPLEX  SMaTRX (2,2) 

0003 

DIMENSION  SCATER11U0,4) 

0004 

COMMON  /rKSCT/SCATER, SMATRX 

OOuS 

GOTO (1 00, 200,300,400) If  IX  ISC ATEN  II# 1>) 

000b 

100 

CALL  PLAIE(I) 

OOoT 

RETURN 

oooa 

200 

CALL  01HEUCIJ 

OOoO 

RETURN 

uoto 

300 

CALL  TKIHEOU) 

001 1 

return 

0012 

400 

CALL  OIPOLE(I) 

0013 

RETURN 

0014 

END 

FORTRAN 

IV 

STORAGE 

MAP 

NAME 

OFFSET 

ATTRIBUTES 

I 

000014 

InTEGER*2  parameter  variable 

IF  I X 

OQuOOO 

1nTEGER*2  PROCEDURE 

plate 

OOOOOO 

R£AL*4 

PROCEOuRE 

uiheu 

UOOOVO 

REAL  a  4 

PRuCEDuRE 

trihed 

OOoOOO 

R£AL*4 

PROCEDURE 

DIPOLE 

OOOOOO 

R£AL*4 

PROCEDURE 

COMMON 

SLOCK  /NKSCT / 

LEnGTh  003140 

SCATER 

OOOOOO 

R£AL«4 

ARRAY  (100,4)  VECTORED 

SmaTkx 

003100 

ComPLEX*S  ARRAY  (2,2)  VECTORED 
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fortran  IV 


V01C-03F *  THU  28-OCT-82  00!07:S1 


0001 

c 

c 

c 

00u2 

ooui 

0004 

OOoS 

oOOb 

0007 

oooa 

0004 

0010 

0011 


subroutine  plate  id 

PLAT  PLATE  SCATTERING  MATRIX 

U I MENS  1 ON  SCATER 1 100 «  4) 

COMPLEX  SM A  T  MX (2/2) 

COMMON  /kyKSCT/SCATER.SMATRX 
SHSIGM=SURT(SCATER(I,2)) 
SmATKX(1, 1)»CMPLX(-1.,0.)*SMSIGM 
SMATRX(1,2)sCMPLX(0.,0.)*SRS1GM 
SmATRX(2, 1)=SMATKX(1,2) 

SMATMX (2,2)=SMATHX(1,  1} 

ME  t  URN 
end 


PAGE  001 


FOkTRAN 

IV 

STORAGE  MAP 

NAME 

OFFSET 

ATT  R IdUTES 

I 

SRSIGM 

SORT 

CMPLX 

000014 

OOwOib 

OOoOOO 

O0O000 

INTEGER*2 
REAL  *  4 

RE  AL*4 

complex**) 

parameter  variable 
variable 

PROCEOURE 

RRUCEOuRE 

COMMON  dLOCK  /«K  SC  T /  LENGTH  00314 0 

SCATtR  000000  REAL*4  ARRAY  (100,4)  VECTORED 

ShaTkX  003100  C0MP0EX*a  ARRAY  (2,2)  VECIUREO 
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■-  v-***#***-*,^  ,  - 


FORTRAN  IV  VO lC«03f ♦  ThU  28-OCT-62  o050?S5l>  P*6E 

C 

oooi  subroutine  oineom 

c 

C  OlHtDRAL  SCATTERING  MATRIX 

c 

U0O2  DIMENSION  scaterciuo,  aj 

0003  COMPLEX  SMATKX12,2) 

OOum  COMMON  /mKSCT/SCATER/SMATRX 

0005  SNSIiiMsSwRI  (SCaTER(1,2)  ) 

0006  SMATRxll, 1)=CNPLX(C0812.»SCATERII, 3) ) ,0.)«SMSIG* 

OOol  SMATRXll #  2) sCMPLX lSINl2,»SCATtRH»3))#o.J*Srt5I&M 

0000  SMATkX (2. 1 ) sSMATMX ( 1 ,  2) 

0009  SMATHX(2,2)sCNPLX(-COSl2.*SCATER(Ir3)),0.)*SKSI&N 

0010  RETURN 

0011  ENO 


fortran 

IV 

STORAGE  HAP 

NAME 

OFFSET 

ATTRIBUTES 

1 

000014 

INTEGEH*2 

PARAMETER  VARIABLE 

SRS16M 

000036 

REAL44 

variable 

SQmT 

OOOOOO 

Rt AL*4 

PROCEDURE 

cmplx 

OOOOOO 

COmPLEX*» 

procedure 

COS 

Oowuuu 

NEAL*4 

PROCEDURE 

SIN 

UOUOOO 

REAL*4 

PROCtDoRE 

COMMON  SLOCK  /HKSCT /  LENGTH  003l«0 

SC A  TER  OOOOOO  REal*4  ARRAv  (100, 4J  VECTORED 

SMAIKX  003100  COmPLEX*S  ARRAY  (2,2)  VECTORED 
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FOHTRAr,  IV 


V01C-03F*  TrlU  26-QCT-82  00S08S02 


C 

OOOl 

C 

c 

c 

uoue 
0003 
uO  o  a 
ouw  5 
OOOb 
OOv7 
uOua 
0009 
0010 
0011 


SUBROUTINE  TRlHEUlI) 

TRIHEDRAL  SCAIIERXK6  MATRIX 

01m£ nSIOn  SC A  TER (100,4) 

COMPLEX  SMaTRX(2,2) 

COMMUN  /«»KSCT/SCATER,SHATRX 
SHSIOmsSwRT (3CATER(I,2) ) 

SmaTrx ( 1 , 1) *CMPLX (-1 .,0.)*SHSIGM 
3MATRX (1,2) *CMPL  X  1 0 . » Q  «  )  *  SRS IbM 
SMA  Tk  * (2, 1 )sSMATmX (1  * 2) 

3M A TRX(2,2)  *3MA  TRX (1,1) 

RETURN 

ENO 


FOR ThAn  IV 


STORAGE  RAP 


name  offset  attributes 


I  O00014 
SRSlfiM  u0o03b 
SORT  UOUOOO 
CMPLX  000000 


INTEGER*2  PARAMETER  VARIABLE 
real**  variable 

real**  procedure 

complex *B  PROCEDURE 


COMMON  BLOCK  /RKSC T /  LE»6TH  003140 

SCAltR  OOuOuO  REAL**  ARRAY  1100,4)  vectored 

SMaTrX  uOSluo  COMPLEX*#  ARRAY  (2,2)  VECTbRED 


FORTRAN  IV 


VO IC"03P ♦  Thu  28-UCT-62  ou:o8!07 


PAGE  00 


C 

0001 

C 

C 

c 

0oo2 

oo03 

0004 

0005 

OOub 

0007 

oooa 

ooov 

0010 

uou 


SUUKUUT  INE  DXPULEU) 

U I POLE  SCATTERING  MATRIX 

DIMENSION  SCATERU00.4) 

COMPLEX  SMATHXC2.2) 

COMMON  /wKSCT/SCaTER.SMATRX 
SRSIUMsSORT  (SCATER  (1,21) 

SMATRx(l.l) sCMPLX C-C0S(2,*SCATER(I»3) ) ,0.) *SRSIGM 
SMATRXU,2)sCMPLX(-COS(SCATERlI,3))*SIN(SCATtR(I,3)),0.)* 

1  SKSIGM 

SMATKx 12# 1 ) sSMATRX ( 1 . 2) 

5MATRX (2,2)sCMPlX(-SIN12,«SCATER  H.3) ) #0.)*SRSIGM 

RETURN 

Enu 


FORTRAN 

IV 

STORAGE 

MAP 

NAME 

OFFSET 

ATTRIBUTES 

I 

000014 

INTEGEMA2 

parameter  variable 

5RS1GM 

000036 

REAL*4 

variable 

SORT 

000000 

RE AL*4 

PROCEDURE 

CMPLX 

OOOOoO 

COMPLEX*# 

PROCEOUHE 

COS 

oooooo 

REAL*4 

PROCEDURE 

SIN 

OGuOuO 

r£Al«4 

PROCEDURE 

COMMON 

block  /rksct/ 

LENGTH  003140 

SCATER 

OOOOOO 

REAL*4 

ARRAY  (100.4)  vECTURED 

SMATRX 

003100 

COMPLEXES 

ARRAY  12.2)  VECTORED 
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FORTRAN  IV 


V01C-03F+  TrtU  28-UCT-82  00:17:32 


PAGE  001 


0001  SUbRUUTlNE  HEAUEK 

C 

c  plotting  heaoer  data  printout 
c 

00o2  COMPLEX  SMATKX(2,2) 

O0u3  DIMENSION  SC  A  TER ( 1 00 , 4) 

0004  INTEGER  IFIL£lb),NXMIT(2),NPNES,AY,ltN 

0005  COMMON  /N0HKF/IFSTFa,IGP,LSTEP,NST£P,0F.CF,F6«» 

OOob  COMMON  /HKSCT/  SC A TER , 5MA TH X 

000/  COMMON  /HEAQ/AISOL, NSC A T, G A INA , NO  I SE , 

1  RANGE, ObLUSS,NXMIT, IF ILE , SO , B I G , 

1  SwRrt,3NKV,SNRHl , SNRHQ , SnR V I , SnR VO , 3NH 
oOoe  COMMON  /SIGNAL/PTPnR,  R ANGE4 , CR ,  ANTG2,SL0SS,P14C 

U0o9  DATA  KY.HN/'  r','  N'/ 

C 

C 

00 10  i»PR£S=KN 

0011  IF ( NO  I SE . EG . 0 ) GO  T  0  5 

0013  NPRESsKY 

C 

0014  5  CALL  PLOTI-l, 0,760) 

0015  CALL  V14CSZI4) 

001  fa  TYPE  <7 

0017  TYPE  10.IF1LE,NSCAT,NSTEP,GAINA,AIS0L,NXMIT,PTP*R, 

1  CR,08L0SS,BIG 

00 Id  TYPE  U.NPRES,  (SO/l.E-b>-,SNRh,SNKV, RANGE 

0019  9  FORMAT (1M>,40X,  'WF  GulOAhCE  TECHNOLOGY  POLARIZATION  SIMULATION') 

0020  10  F  OHM A  T (9X# 'OATA  FILE  NAME  S ' • 6A2  # 

1  lX.'NUM.  SCATTEkEKS: *, 13, 

1  3X.'FREG  STEPS: ' , 13, 

1  JX ,  ' AN  T  GAIN (Ob ) : ' ,F7,2, 

1  3X,  ' ANT  I30LATI0NI08) ',F7.2Z 

1  9X, 'XMITS ',2A2, 

1  3 X , '  X M I T  Prr/ChnL(*ATTS) : 'F&.2, 

1  3X , 'COMP  RATIO: ' ,Ffa.2, 

1  SX, 'SYSTEM  LOSSIOb) :’ ,F5, 2, 

1  3X , ' FF  T  SCALER:', 1PE14. fa) 

0021  11  format (9X, 'NOISE: ', A2, 

l  3x, 'NOISE  S0IUV0LTS):»,F7.5, 

1  3x,'n  AVG  SnR 108) S ' »Ffa«2» 

1  2X, • V  AVG  SNH(Ob) : ' ,Fb.2, 

1  3X, 'RANGE  TO  TARGET  CELL (ME TERS) : ' , F « .2 ) 

0022  00  20  121,1000 

0023  20  CONTINUE 

0024  CALL  V14CSZI1) 

0025  RETURN 

002b  ENU 
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FORTRAN 

IV 

STORAGE  MAP 

name 

OFFSET 

AT  TRIBUTES 

NPkES 

00064b 

INTEGEK*2 

VARIABLE 

KY 

000014 

INTEGER«2 

VARIABLE 

KM 

000016 

INTEGERS 

variable 

PLuT 

000000 

REAL*4 

PROCEDURE 

v 14CS2 

OOVOOO 

KEAL»4 

PROCEDURE 

I 

000650 

INT£G£R*2 

VARIABLE 

COMMON  dLOCK  /nORKF /  LENGTH  0UQ024 


1FSTFQ 

000000 

IN  T  ESER»2 

VARIABLE 

I  UP 

000002 

INTEGER*2 

VARIABLE 

LSTEP 

000004 

INTEGER»2 

variable 

nstep 

000006 

INTEGER*2 

variable 

OF 

000010 

HEAL*4 

variable 

CF 

000014 

RE  AL  *  4 

variable 

Fh* 

000020 

kEAL*4 

variable 

common 

BLOCK  /I 

RKSC  (/ 

LENGTH  Oo 3 1 40 

SCATtR 

OOOOoO 

REAL*4 

ARRAY  1100,4)  VECTOR! 

SMATRX 

003100 

COMPLEX  *B 

ARRAY  (2,2)  VECTORED 

COMMON 

BLOCK  /HEAD/ 

LENGTH  000114 

A  I  SOL 

OOOOOO 

REAL»4 

VARIABLE 

NSC  A  T 

000004 

INTEGER*2 

variable 

gaina 

000006 

HE AL*4 

variable 

,uISt 

O0O0  12 

integer**! 

variable 

K  A„Gt 

000014 

REAL*4 

V  Ah  I  ABLE 

UHlOSS 

000020 

RE  AL  *  4 

variable 

UmIT 

0  0  u  024 

integer^ 

ARRAY  (2) 

IFILt 

oOoOiO 

1nT£GER*2 

ARRAY  (B) 

SO 

000050 

«EAL«4 

VARIABLE 

BIG 

oooo^a 

RE AL*4 

VAkIAHLE 

SNRM 

OOOOoO 

REAL*4 

VAR  1  ABlE 

SNHV 

000064 

RE AL*4 

VARIABLE 

SNkm  1 

000070 

R£AL»4 

variable 

snrhu 

000074 

hEAL*4 

VARIABLE 

SNRVI 

oooluo 

REAL*4 

variable 

SNkVU 

000104 

REAL*4 

variable 

SNh 

000110 

REAL*4 

variable 

t'lr'-UN 

BLOCK  /SIGNAL/ 

LENGIH  OoOO 30 

P  TPftk 

0  0  0  0  0  0 

R  £  AL  *  4 

variable 

KAMGE4 

0000O4 

REAL*4 

VAR  1  able 

C* 

0  0  0  0  1 0 

RtAL*4 

variable 

AN  TG2 

000014 

REAL*4 

variable 

SLuSS 

000020 

RE  AL  *4 

VARIABLE 

PI4C 

o0o024 

REAL*4 

variable 
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